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Caminante, son tus huellas el camino, y nada mas; 
Caminante, no hay camino, se hace camino al andar. 
Al andar se hace camino y al volver la vista atrds 
se ve la senda que nunca se ha de volver de pisar. 
Caminante, no hay camino sino estelas en la mar. 

Antonio Machado 

1. Introduction 

Natural processes proceed simultaneously on many different scales interacting with 
each other. This is especially important in astrophysics, where an astronomical ob- 
servation or a hypothesized system can only be understood by modeling processes 
on scales spanning many orders of magnitude. Thus, also nuclear physics and as- 
trophysics are closely entwined. Nuclear reactions power quiescent burning of stars 
and cause the most powerful explosions known. They not only release or transform 
energy but also change the composition of the matter in which they occur and thus 
are responsible for the range of chemical elements found on our planet and through- 
out the Universe. Finally, to understand the properties of matter at extreme density 
and/or density, the nuclear equation of state has to be known. It determines the 
properties of neutron stars and the latest stages of the life of stars with more than 
S>Mq which end their lives in a supernova explosion. 

It is evident that nuclear physics input is essential for many astrophysical mod- 
els and this fact is represented in the field of Nuclear Astrophysics. There are 
different perceptions on how to define this field. Some limit it to the application of 
nuclear physics to reactions of astrophysical interest. A more comprehensive, and 
perhaps more adequate, definition includes the more astrophysical aspects in the 
investigation of nuclear processes and nucleosynthesis in astrophysical sites through 
reaction networks. In any case, the interests of astrophysics emphasize different as- 
pects than those of basic nuclear physics and this makes Nuclear Astrophysics a 
distinct research area. From the nuclear point of view, astrophysics involves low- 
energy reactions with light projectiles on light, intermediate, and heavy nuclei. 
Although nuclear physics has moved to higher energies in the last decades, low- 
energy reactions are not well enough explored and still offer considerable challenges 
to both experiment and theory, even for stable target nuclei. Explosive conditions 
in astrophysics favor the production of nuclei far off stability and reaction rates 
for these have to be predicted across the nuclear chart. This proves very difficult, 
especially for first-principle methods, due to the complexity of the nuclear many- 
body problem. The calculation of astrophysical reaction rates also includes special 
requirements and processes not studied in nuclear physics so far. Among these are 
effects appearing in low energy, subCoulomb reactions and reactions on excited 
states of target nuclei. Depending on the conditions, plasma effects also have to 
be considered because these alter the reaction rates. This includes, e.g., the shield- 
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ing of Coulomb barriers through electrons, pycnonuclear burning in the lattice of 
a high-density plasma, and the modification of nuclear partition functions at very 
high plasma temperatures. 

This paper focuses on nuclear reactions in astrophysics and, more specifically, 
on reactions with light ions (nucleons and a particles) proceeding via the strong 
interaction. It is intended to present the basic definitions essential for studies in 
nuclear astrophysics, to point out the differences between nuclear reactions taking 
place in stars and in a terrestrial laboratory, and to illustrate some of the challenges 
to be faced in theoretical and experimental studies of those reactions. The sensi- 
tivity of the reaction rates to the uncertainties in the prediction of various nuclear 
properties is explored and some guidelines for experimentalists are also provided. 
The discussion revolves around the relevant quantities for astrophysics, which are 
the astrophysical reaction rates. The impact of using different models or data is 
always presented with respect to the possible modification of the reaction rates, 
not the cross sections. 

At first, the basic equations through which nuclear processes enter astrophysical 
models are introduced in Sec. [21 The astrophysically relevant energy ranges are 
defined in Sec. [H The thorough discussion of the special stellar effects affecting 
reactions in Sec. |4] is the heart of this review. Stellar and effective cross sections 
are derived in Sec. 14.11 then the relation between laboratory reactions and those 
occurring in a stellar plasma is investigated in Sec. 14.31 As it is important to realize 
when individual reactions are important and when they are not, reaction equlibria 
are introduced in Sec. 14.21 Cross sections in a plasma are affected by the free 
electrons present which shield the nuclear charge. This is explained in Sec. 14.41 
Finally, astrophysically relevant reaction mechanisms are reviewed in Sec. [S] and 
reactions through isolated resonances (Sec. 15. 2| ). in systems with high nuclear level 
density (Sec. I5.3.2p . and direct reactions (Sec. 15. 5p are discussed separately. This 
includes a detailed discussion of the sensitivities of the rates on the required input 
in Sees. I5.4T2] and [5.5.41 Section [6| presents a brief conclusion. 

2. Reaction networks and rates 

Nuclear reactions are the engine of stellar evolution and determine the overall pro- 
duction of the known chemical elements and their isotopes in a variety of nucle- 
osynthesis patterns. A detailed understanding of the characteristic production and 
depletion rates of nuclei within the framework of a nucleosynthesis process is crucial 
for reliable model predictions and the interpretation of the observed abundances. 

Instead of the number AI of nuclei of a given species per volume V (the number 
density n = M/V), it is advantageous to use a quantity independent of density 
changes: the abundance Y = n/(ppiaiVA), where ppia is the plasma density and 
TVa denotes Avogadro's number. The change of abundances Y with time due to 
nuclear processes is traced by coupled differential equations. Due to the nature of 
the involved reactions and the vastly different timescales appearing, the equation 
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system is non-linear and stiff. In addition, for complete solubility of the coupled 
equations the number of equations N has to equal the number of involved nuclei 
acting as reaction partners and thus an equation matrix of size has to be solved. 
Nucleosynthesis processes include thousands of nuclides and tens of thousands re- 
actions. This still makes it impossible to fully couple such a reaction network to a 
full set of hydrodynamic equations as would be required for a complete modeling 
of nucleosynthesis in a given astrophysical site. 

A reaction network generally can be written as ^ 

PplaiVA PplaA'A ^ ^ ^ J 

(1) 

where 1 < i < numbers the nucleus, iXj is the jth rate for destruction or creation 
of the ith nucleus without a nuclear projectile involved (this includes spontaneous 
decay, lepton capture, photodisintegration), and iVj is the rate of the jth reaction 
involving a nuclear projectile and creating or destroying nucleus i. Similarly, we have 
three-body reactions where nucleus i is produced or destroyed together with two 
other (or similar) nuclei. Reactions with more participants (denoted by . . . above) 
are unlikely to occur at astrophysical conditions and are usually neglected. The 
quantities jKj, ^Kj, and f A'j/c are positive or negative integer numbers specifying 
the amount of nuclei i produced or destroyed, respectively, in the given process. As 
shown below, the rates A, r, and f contain the abundances of the interacting nuclei. 
Rates of type A depend on one abundance (or number density), rates r depend on 
the abundances of two species, and rates r on three. 

The rates iXj appearing in the first term of Eq. ^ are reactions per time and 
volume, and only contain the abundance Yj. For example, iXj is simply rijLj — 
YjPpiaNALj for /3-decays. The factor Lj = (ln2)/-'Ti/2 is the usual decay constant 
(with the unit 1/time) and is related to the half-life -'T1/2 of the decaying nucleus 
j. It has to be noted that some decays depend on the plasma temperature and thus 
Lj is not always constant, even for decays (see Eq. (|29|) in Sec. 14. ip . 

Two-body rates r include the abundances of two interacting particles or nuclei. 
In general, target A and projectile a follow specific thermal momentum distributions 
duA and dua in an astrophysical plasma. With the resulting relative velocities 
va — Va, the number of reactions per volume and time is given by 

rAa = J a{\vA - Va\)\vA - ValduAdUa , (2) 

and involves the reaction cross section tr as a function of velocity, the relative 
velocity va — Va and the thermodynamic distributions of target and projectile duA 
and dUa- The evaluation of this integral depends on the type of particles (fermions, 
bosons) and distributions which are involved. 

However, many two-body reactions can be simplified and effectively expressed 
similarly to one-body reactions, only depending on one abundance (or number 
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density) . If reaction partner a is a photon, the relative velocity is always c and the 
quantities in the integral do not depend on dua- This simplifies the rate expression 
to 



where Lj{T) stems from an integration over a Planck distribution for photons of 
temperature T. This is similar to the decay rates introduced earlier and therefore 
we replaced r by A in our notation and can include this type of reaction in the first 
term of Eq. ((T)) . A similar procedure is used for electron captures by protons and 
nuclei. Because the electron is about 2000 times less massive than a nucleon, the 
velocity of the nucleus is negligible in the center-of-mass system in comparison to the 
electron velocity (|wnucicus— ^Vioctronl ~ | ■^'electron | ) • The electron capture cross section 
has to be integrated over a Fermi distribution of electrons. The electron capture 
rates are a function of the plasma temperature T and the electron number density 
Ue = YeppuNA- In a neutral, completely ionized plasma, the electron abundance 
Ye is equal to the total proton abundance Ye = J2i ^i^i ^^'^ thus 



Again, we have effectively a rate per target L (with unit 1/time) similar to the 
treatment of decays earlier and a rate per volume including the number density 
of only one nucleus. We denote the latter by A and use it in the first term of 
Eq. ([l|). This treatment can be applied also to the capture of positrons, being 
in thermal equilibrium with photons, electrons, and nuclei. Furthermore, at high 
densities (ppia > lO^^gcm^'') the size of the neutrino scattering cross section on 
nucleons, nuclei, and electrons ensures that enough scattering events occur to lead to 
a continuous neutrino energy distribution. Then also the inverse process to electron 
capture (neutrino capture) can occur as well as other processes like, e.g., inelastic 
scattering, leaving a nucleus in an excited state which can emit nucleons and a 
particles. Such reactions can be expressed similarly to photon and electron captures, 
integrating over the corresponding neutrino distribution. 

In the following, we focus on the case of two interacting nuclei or nucleons as 
these reactions will be extensively discussed in the following sections. (We mention 
in passing that Eq. ([2]) can be generalized to three and more interacting nuclear 
species by integrating over the appropriate number of distributions, leading to rates 
f and higher order terms in Eq. p]l.) 

The velocity distributions in the rate definition in Eq. ([2]) can be replaced by 
energy distributions. Furthermore, it can be shown th at th e two distributions can 
be replaced by a single one in the center-of-mass system.EEl Then the two-body rate 
r is defined as an interaction of two reaction partners with an energy distribution 
(j){E) according to the plasma temperature T and a reaction cross section cr'(E) 



Xa = Lj{T)nA 



(3) 



A, 



'nucleus, CC ~ 




(4) 
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specifying the probability of the reaction in the plasma: 



rAa = ^^ a'^SE)4>{E)dE . (5) 

1 + OAa J 


The factor 1/(1 + Sau) with the Kronecker symbol S is introduced to avoid double 
counting. The nuclear cross section aAa is defined as in standard scattering theory 

by 

number of reactions target~^sec~^ , , 

(TAa = 5 -r- ■■ ^—r, ■ (6) 

nux of incoming projectiles 

However, in an astrophysical plasma reactions not only proceed on the ground state 
of a nucleus but also from excited states. This is implied in the notation for the 
modified cross section ct^^, contrasting the usual laboratory cross section (denoted 
by (TAa or cr^Aa ) reactions acting only on the ground state of the target nucleus. 
The implications of using such a cross section modified in the stellar plasma, instead 
of the usual laboratory one, will be discussed in Sec. 14.11 It should be noted that 
a' may be a function not only of energy but also of plasma temperature. 

The distribution of kinetic energies of nuclei in an astrophysical plasma with 
temperature T follows a Maxwell-Boltzmann distribution (MBD) (/){£) = 0mb(T) 
and we obtain finally: 

"■^^ = ^ r4F.P^,,^^l(a^;)^, (7) 

1/2 °° 

{<^v)Aa=(^) {kT)-^'^ I Ea'^,{E)e-^ dE . (8) 
yrriAaTr J J 


Here, toao denotes the reduced mass of the two-particle system and {<yv)Aa 
action rate per particle pair or reactivity under stellar conditions. The angle brackets 
stand for the appropriate averaging, i.e. integration, over the energy distribution. 
For the remainder of the paper we will be concerned with the determination of this 
reactivity and the involved cross sections, respectively. 



3. Relevant energies 

Before we proceed to the details of the determination of the reaction cross sections, 
it is instructive to further investigate the rate equation and to derive the relevant 
energies at which the nuclear reaction cross sections have to be known. Although 
the integral in Eq. ([8]) runs to infinity, the MBD folded with the cross section selects 
a comparatively narrow energy range with non-negligible contributions to the total 
value of the integral. Historically, this energy range is called the Gamow window 
because Gamow realized early on the astrophysical relevance of the fact that - if 
the energy dependence of the cross section is dominated by the Coulomb barrier 
between the projectile and the target - the integrand in Eq. ([8|) can be factorized 
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ax 

F = Ea{E)e-^ = S{E)e--^e^^ , (9) 

where S is the astrophysical S'-factor 

S = aEe^ (10) 

which is assumed to be only weakly dependent on the energy E for non-resonant 
reactions. The second exponential in Eq. (|9]) is called the Gamow factor and contains 
an approximation of the Coulomb penetration with the Sommerfeld parameter 



' h y 2E ' ^ ' 

where Za, Za are the charges of projectile a and target A, respectively, and ruAa 
is their reduced mass. While the first exponential (the tail of the MBD) decreases 
with increasing energy, the Gamow factor increases, leading to a confined peak of 
the integrand, the so-called Gamow peak. The location of the peak Eq is shifted 
to higher energies with respect to the maximum of the MBD at Emb — kT {kT — 
Tg/l 1.6045 MeV when Tg is the plasma temperature in GK). The width of the peak 
gives the astrophysically relevant energy range in which most of the reactions will 
take place. 

In absence of a Coulomb barrier the energy dependence of the non- resonant cross 
section is roughly given by the one of the wave number of the particle (1/ Ve) folded 
with the angular momentum barrier. This does not, however, lead to a relevant shift 
of the peak of the integrand compared to the peak of the MBD. Thus, the effective 
energy window for neutrons is simply the peak of the MBD. 

The above considerations concerning the location and size of the energy window 
have given rise to simple approximation formulae extensively used by experimen- 
talists to estimate the energies of interest. For instance, with a charged projectile 
the location Eq and width A of the Gamow window is often computed from 

Eo = 0.12204 (mAaZ^aZATg) ' , (12) 

1 

A = 0.23682 (mAaZlZlTi) ' , (13) 

which is derived from Eq. ([9]) assuming a Gaussian shape of the peak and using 
appropriate numerical constants, yielding Eq and A in MeV.^^EIll ^j-^jg approxima- 
tion of the Gamow window Eq ± (A/2) is valid for some but not all cases because 
it is oversimplified. The above factorization with the given Sommerfeld parameter 
including the charges of projectile a and target A implicitly assumes that the energy 
dependence of the cross section is given by the Coulomb penetration in the entrance 
channel a + A of the reaction a + A ^ B + b. It has been realized, however, that 
sometimes resonances below the Gamow window derived with the above appr oxi- 
mation significantly contribute to the reaction rate for certain capture reactions .^^^^ 
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integrand 
standard formula 




E [MeV] 



Fig. 1. Comparison of the actual reaction rate integrand and the Gaussian approximation of 
the Gamow window for the reaction ^^^Sn{p,a)^^^ln at T = 5 GK. The two curves have been 
arbitrarily scaled to yield similar maximal values. The maximum of the integrand is shifted by 
several MeV to energies higher than the maximum Eq of the Gaussian. (Taken from Ref. 6 with 
kind permission.) 



In those considered cases the energy dependence of the cross section is dominated 
by the energy dependence of the 7 width in the exit channel instead of the charged 
particle width in the entrance channel. This can be generalizecP and leads to the 
important realization that the energy dependence of the integrand T has to be 
numerically examined in order to derive reliable energy windows. This has been 
generally and extensively studied in Ref. [5] Here, only a few examples are shown. 
Figure [1] shows a comparison between the actual T and the integrand assumed with 
the standard approximation. In this case, the relative shift of the energy window is 
to higher energy because of the higher Coulomb barrier in the exit channel. In other 
cases, the shifts can also be to much lower energy than predicted by the standard 
approximation. 

For a detailed understanding of these differences one has to realize that both 
resonant and Hauser-Feshbach cross sections (see Sees. 15.21 15.3.21) can be expressed 
aJ3 

a(x^(2J„ + l)^i^ , (14) 

with X being either Breit-Wigner widths or averaged Hauser-Feshbach widths, 
depending on the context. The width of the entrance channel is given by , 
the one of the exit channel by X^" , and the total width including all possible 
emission channels from a given resonance or compound state with spin J„ by X^^l — 

xi" + xi" + ... 

It has become common knowledge that a cross section of the form shown in 
Eq. (|14p is determined by the properties of the smaller width in the numerator if 
no other channels than the entrance and exit channel contribute significantly to 
X^q'J (see also Sec. I5.4.T|1 . Then X^^^ cancels with the larger width in the numer- 
ator and the smaller width remains. (The effect is less pronounced and requires 
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a more detailed investigation when other channels are non- negligible in X^^^.) In 
consequence, the energy- dependence of the cross section will then be governed by 
the energy dependence of this smallest X'^ . Only if this happens to be the charged- 
particle (averaged) width in the entrance channel, the use of the standard formula 
for the Gamow window (Eqs. [T^ fT5|) will be justified. Since X-^^ and X^ have dif- 
ferent energy dependences, it will depend on the specific energy (weighted by the 
MBD) which of the widths is smaller. Therefore, the above approximation should 
not be applied blindly but rather the actual Gamow windows have to be deter- 
mined from the true energy dependence by inspection of the integrand in Eq. ([8]). 
Extensive tables of revised effective energy ranges for astrophysics from such a 
numerical inspection are given in Ref. [6] These ranges can be shifted by several 
MeV to higher or lower energy compared to the ones obtained with the standard 
formula. Furthermore it is found that the assumption of a Gaussian shape of IF 
is untenable for the majority of cases with intermediate and heavy target nuclei. 
Rather, the integrand !F may show a pronounced asymmetry around its maximum 
value. Therefore, the energy of the maximum alone is not sufficient to determine 
the astrophysically relevant energy range. 

Although derived from cross sections of a specific model prediction, the energy 
windows given in Ref. [6| are supposed to be robust. This can be understood by 
realizing that they mainly depend on the relative energy dependence of the acting 
reaction channels and not the absolute value of the cross sections. This dependence 
is governed by the relative energy and the Coulomb barrier seen in each reaction 
channel. Therefore, the limits of the energy windows are set by the knowledge of the 
charges of the nuclei involved in the different reaction channels, and the reaction 
Q- values. Only the latter may be unknown for nuclei far off stability and mass 
measurements may have an impact. 

Another important consequence of using the correct energy dependence is that 
different reactions may not have necessarily the same effective energy window, even 
when projectile and target nucleus are the same. This is immediately seen when 
considering the case of a reaction with a positive Q-value for one channel but a 
negative one for another reaction channel. The astrophysically relevant energy win- 
dow of the exothermic reaction may lie below the threshold of the other reaction 
channel. Obviously, the relevant energy window for the endothermic channel cannot 
open below the threshold energy and thus has to lie at higher energies than the 
one for the capture. A randomly chosen example for such a reaction pair would be 
i°''Pd(p,7)^°^Ag and i°'^Pd(p,n)i°4Ag^ with the energy windows at 1.5 - 2.8 MeV 
and 5.07 — 5.7 MeV, respectively, for a plasma temperature of 2 GK.^ Below the 
(p,n) threshold at 5.06 MeV, the (p,n) cross section is zero and does not give a 
contribution to the integral in Eq. ([8|). This also further illustrates the limitation 
of the standard approximation which yields identical energy windows for reactions 
with identical entrance channels but different exit channels. It has to be noted, 
though, that the effective energy windows only point out the energy range con- 
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tributing mostly to the rate integral at given stellar temperature but do not make 
a statement on the size of the rate or its astrophysical relevance. 

It has already been mentioned that Eq. (|T4)) applies to reactions either exhibiting 
isolated resonances treatable by a Breit-Wigner resonance formula or smooth cross 
sections stemming from an averaging over a large number of narrowly spaced and 
unresolved resonances. Therefore the derived energy windows are also applicable to 
obtain the relevant energy ranges in which narrow resonances have to be considered. 
They do not, however, specify the relative strengths of the resonances within a given 
window. 

Definition (fTO|) only makes sense when using the laboratory cross section cr = 
^lab ^ggg ([T51) in the following section). Strictly speaking, with the laboratory 
cross section the energy windows apply to laboratory measurements only, i.e. to 
the determination of the ground state component of the actual stellar cross section. 
For low stellar temperatures and positive Q-values (see the discussion in Sec. 14. 1|) . 
however, this will dominate the stellar cross section. Relevant energy windows can 
also be derived numerically for stellar cross sections (see next section), of course, 
in the same manner. Since these cannot be measured (yet), they may be of limited 
use, though. 

A general scrutiny of the astrophysically relevant energy windows up to a plasma 
temperature of 5 GK (above this temperature, reaction equilibria are established 
which do not require the knowledge of individual rates; see, e.g., Refs. [TJ ^ shows 
that the appearing interaction energies are small by nuclear physics standards. For 
neutron-induced reactions, the encountered maximum energies are a few hundred 
keV, depending on the e xamin ed nucleosynthesis process (e.g., in the s process they 
are more like 8-60 keV) .UJ^EI The relevant energy windows are shifted to higher 
energies for charged reactants, with a few MeV for reactions with protons and 
several MeV up to about 10 MeV for reactions involving a particles. The formulae 
given in Eqs. p2)) . (fT3|) are inadequate for the determination of these energies and 
should not be used anymore. 

4. Stellar effects 

4.1. Stellar cross sections and reciprocity of stellar rates 

In an astrophysical plasma, nuclei quickly (on the timescale of nuclear reactions 
and scattering) reach thermal equilibrium with all plasma components. This allows 
thermal excitation of nuclei which follows a Boltzmann law and gives rise to the 
stellar cross section 



where the sum runs over all excited states ^ of the target nucleus A (for simplicity, 
here we assume the projectile a, i.e. the second reaction partner, does not have 




(15) 
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excited states) with spin factor = 2J^ + 1 and excitation energy E^. Thus, the 
stellar cross section is the sum of cross sections (evaluated at their respective 
center-of-mass energies ^ E — E^ with = when Ef^ < 0) for reactions 
on a nucleus in excited state fi, weighted by the population coefficient 7-"^. The 
above relation can be derived from a Saha equatiorPl. It is to be noted that the 
stellar cross section depends on energy and temperature, contrary to the usual cross 
section which only is a function of energy. 

The quantity G is the partition function of the nucleus. Often, the partition 
function normalized to the ground state 



Go(T) - 




is used (the ground state is labeled as /i = 0, the first excited state as fj, = 1, . . .). 
Equation p6p shows how the computation can be extended beyond the energy of the 
highest known discrete energy level ^^^^^ by using an integration of a nuclear level 
density p over a range of excitation energies e. Likewise, the sums appearing in (jl5l) 
can be amended with a supplemental integration over the level density above the last 
discrete level used. Although the product of the Boltzmann factor (2J+ l)e~'^/*^'^^^ 
and p does not have a trivial energy dependence, it has been shown that for the 
application of (|T6)) at temp eratures T < 10 GK it is sufficient to integrate only 
up to iJmax = 25 MeV.l211QI Temperatures above 10 GK are encountered in some 
explosive astrophysical events, in accretion disks, and in the formation of neutron 
stars and black holes. Nuclear transformations in such environments are described 
in reaction equilibria between several or all possible reactions (with the exception 
of reactions mediated via the weak interaction because they are too slow in most 
cases), replacing full reaction networks by simplified abundance equations (see, 
e.g., Refs. [U |3]for details) but still containing the partition functions. At such high 
temperatures, a straightforward application of ()16|) would overestimate the partition 
function because continuum effects have to be taken into account. These can be 
treated by appr oxim ated correction factors to p and extending the integration to 
Eiaax ^ 35 MeV.fl^ A more rigorous treatment of the correction would be desireable 
but the sheer number of involved nuclei proves prohibitive for fully microscopic 
approaches. 

The use of stellar cross sections in the calculation of reaction rates assures an 
important property, the reciprocity of forward and reverse rate. A scheme of the 
energetics and the transitions between nuclear levels in the involved nuclei is shown 
in Fig. [2] For a reaction a + A ^ b + F involving only one initial level p in nucleus 
A and one final level ly in nucleus F (this can also be the compound nucleus C 
if the ejectile & is a photon) the well-known reciprocity relation between forward 
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Fig. 2. Schematic view of the transitions (full arrows denote particle transitions, dashed arrows 
are 7 transitions) in a compound reaction involving the nuclei A and F, and proceeding via a 
compound state (horizontal dashed line) with spin Jq and parity tTq in the compound nucleus C. 
The reaction Q values for the capture reaction (Qcap) and the reaction A— >F (QF=QAa) are given 
by the mass differences of the involved nuclei. Above the last state, transitions can be computed 
by integrating over nuclear level densities (shaded areas). 



i'^Aa) ^^"^ reverse (cr^^) reaction cross section iJSHI] 

'^FbiEF) = -y —a'XaiEA) , (17 

where rriyi, mp are the reduced masses, i?^, Ep the center-of-mass energies relative 
to the levels /i, z/, respectively, and g the spin factors as before. This relation 
connects one initial with one final state and therefore is not applicable to the 
regular laboratory reaction cross sections which connect one initial state (usually 
the ground state /i = 0) of the target nucleus with a number of possible final states 

a = a'^^=a^=' = a'=Y,aX ■ (18) 

In order to obtain a cross section obeying reciprocity one has to construct a theoret- 
ical quantity called effective cross section a^^ . The effective cross section is a sum 
over all energetically possible transitions of initial levels to final levels (capture: in 
nuclei A and C, otherwise: in nuclei A and F; as indicated in Fig. [2]), applying to 
compound reactions as well as to direct reactions (see Sec. [5]). It includes all the 
transitions shown by arrows in Fig. [2] and therefore sums over all final levels v and 
initial levels /i 

<(i?o) = EE-#'-"^(^.)-E-§'-"(^.) • (19) 

90E0 ^ go Eo 
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As before, the relative center-of-mass energy of a transition proceeding from level 
jjL is denoted by E^^ ^ Eq — E^, and cj^" = cr^ = for < 0. The first summand 
(/Lt = 0) in the sum over fi is just the laboratory cross section cr°. Note that the 
effective cross section is only a function of energy, like the usually defined cross 
section, and does not depend on temperature, contrary to the stellar cross section. 

When interchanging the labels fi and v, a similar quantity ct^ is obtained for 
the reverse direction commencing on levels v. It is straightforward to show that the 
two effective cross sections obey the reciprocity relations 

A rpA 

off 9o 9a mAaJ^o off /r,^^ 

95 9b mpbES 

A 

„eS rpF 9o 9a rriAa eff j^A /oi \ 

(^FbEo - — (^Aa^o , 

95 9b mFb 



which are identical to the one in ([T7| for a single transition between two states in 
two nuclei. The relative energies of the transitions proceeding on the ground states 
of the two target nuclei for forward and reverse reaction are denoted by Eq and 
Eq, respectively. 

The effective cross section is, of course, unmeasureable. Its usefulness becomes 
apparent when we combine definition ()15|) of the stellar cross section with the 
definition (jS]) for the astrophysical reaction rate. Since excited states of target nuclei 
are populated in a stellar plasma according to ()15p . we have to sum over the rates 
for reactions from each level and weight each summand with the population factor 



i— 



1/2 



-3/2 



dEf. 



(22) 



This means that projectiles with MB distributed energies are acting on each level 
/Lt separately. Insertion of definition for the population factor leads to 



ioo 1 oo 

/ <aE^^-^ dE^ \=T.[ VAaE^e^^ dE^ 
J ^ 



M 



--EJ/(fcT) 1 p A E'^ + E'^ 

^-^ a^.E^e-^ dE^ = \^ ilL^^^E^e-^ dE^ 



M 

(23) 



where Gq is the normalized partition function as defined in (|16p . In order to obtain 
an expression similar to the original single MBD, the integral can be transformed 



June 10, 2011 0:11 WSPC/INSTRUCTION FILE rauscher 



The Path to Improved Reaction Rates for Astrophysics 15 

by replacing dE^ — dE^ , with E^ = E^ + E^ and this yields 



oo ^ 

= ^ E / '^<aiE^ - - KV-^ dE^ = 

OO ^ 

^^aT. f'^<aiE^)[E^-K]^-^dE^ ■ (24) 







In the last line above, the lower limit of the integration was reset to Zero. This is 
allowed because cross sections at negative energies do not give any contribution to 
the integral. It has been pointed out in Ref . fT2l t hat it i s mathematically equivalent 
when sum and integral are exchanged, leading tj ^^ l l^ l 



, 1/2 

!Aa ^ 



-) m-^'^ E / K - K] dE^ = 



The last line was obtained by realizing that the expression in the curly brackets is 
identical to o'^^^E^ , with the effective cross section from ([19]). Thus, the weighted 
sum over many MBDs acting on the thermally populated excited states is reduced 
to a single MBD acting on an effective cross section and divided by the normalized 
partition function. In terms of relevant physics, this means that the Boltzmann 
factor in the population probability is offset by shifting down each MBD to the 
same relative energy.tlSI Now a' in ([8|) can be identified as a' = cr°^ /Gq, introducing 
a temperature dependence while cr°^ is conveniently independent of T. 

It may be confusing that the stellar reactivity frequently is written as (cr*w) 
instead of (crv)* . This is not meant to imply that the stellar cross section as defined 
in p5|l is inserted in a single integral over a MBD as shown in ([8|). Rather, the 
angle brackets imply a separate integration for each populated state as performed 
in ([22]) and ([23l) in this case. 

Equation (|25p not only simplifies the numerical calculation of the stellar rate but 
also allows to better understand certain details. For instance, it can immediately 
be seen that stellar rates obey a reciprocity relation because the effective cross 
sections do. A similar expression has to hold for the reverse reactivity {crv)*pf^ as for 
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the forward reactivity {<yv)\^, being 

1/2 



Mk=f^) (fcr)-3/^-L/afXe-*d< . (26) 
ympbi^ J Go J 



This expression is derived in the same manner as psp but by starting from thermally 
populated excited states in the final nucleus. With the help of (1^ w e can express 
the reverse reactivity of (|26|) in terms of the forward reactivityl^^^^l 



)*Fb ^ 9n_9a Gq ^ niAa ^ p-QAa/(kT) 
)*Aa 



9o9b \mBbJ 



where QAa = Eq — Eq is the reaction Q-value of the forward reaction. The reci- 
procity relation (j2ip applies to photodisintegration and captures as well. In re- 
lating the photodisintegration rate A to the capture rate, however, it has to be 
assumed that the denominator exp(£^/ (kT)) — 1 of the Planck distribution for pho- 
tons appearing in the photodisintegration reactivity L-y — {<Jv)'^,^ can be replaced 
by exp{E/{kT)), similar to the one of a MBD with the same temperature T. With 
this approximation and realizing that g~f = 2, one obtaina^iSilSiMl 



(28) 



Ej _ 9o9a Gq f mAakT \^^ -QAa/ikT) 

{'^v)*Aa " 9S G§ { 2nh^ J ' 

in the same manner as (j27p . Using the approximation of the denominator resulting 
in (|28p is very important for the application in reaction networks. Employing the 
expressions (|26p and (|28p avoids numerical inconsistencies in network calculations 
which may arise when forward and reverse rates are calculated separately (or even 
from different sources) . The proper balance between the two reaction directions can 
only be achieved in such a treatment. Furthermore, simplified equations for reaction 
equilibria (see Sec. 14.21) can be derived which prove important in the modeling and 
understanding of nucleosynthesis at high temperature. 

How large is the error stemming from the approximation involved in the 
derivation of ([28| ? Although mathematically unsound, it turns out that setting 
exp{E / {kT)) — 1 w exp{E / (kT)) is a good approximation for the calculation of the 
rate integrals and introduces an error of less than a few percent for astrophysically 
relevant temperatures and rate values .^^CSllllS] jjj other words, the contributions 
to the integral in ([3]) are negligible at the low energies where Planck and MBDs 
differ considerably (see also Fig. 3.5 in Ref. |3l and Ref. [TH [15]). This is assured by 
either a sufficiently large and positive QAa, which causes the integration over the 
Planck distribution to start not at Zero energy but rather at a sufficiently large 
threshold energy, or by vanishing effective cross sections at low energy due to, e.g., 
a Coulomb barrier. The assumption may not be valid for s-wave neutron captures 
with very small (of the order of Q < kT) or negative Q-valiies, but the required 
correction still is only a few % as can be shown in numerical comparisons between 
photodisintegration rates calculated with the two versions of the denominator. Such 
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Fig. 3. Target nuclei for which an error of 5—10% is introduced in the relation (p,7)/(7,p) at T = 1 
GK with the standard approximation of the Planck distribution. 

a comparison was performed with the code SMARAGD (version 0.8s; see Sec. 15.4.51) 
and the results are shown in Figs. [31 131 and [3 Generally, larger errors appear at 
lower temperature. This results in astrophysical irrelevance of the errors in many 
cases because either the rates are too slow (especially for rates involving charged 
projectiles) or the target nuclei in question are so short-lived that they will never 
be produced at low plasma temperature. The largest error found was between 50 
and 100% for a few heavy nuclei at the driplines for proton- or a-capture at T < 0.3 
GK. For neutron captures, the errors when applying the standard approximation 
for the reverse rate were never larger than 10% at any investigated temperature, 
even at the driplines. In the figures, errors of 5% and smaller are not presented 
in detail because they are assumed to be negligible, especially given the remaining 
uncertainties in the prediction of the rates far from stability. 

Figure |3l shows target nuclei for proton capture where errors reach 5-10% at 
T = 1 GK when computing the photodisintegration rate from ([SS)) . As expected, 
this occurs close to the dripline where the reaction Q- value is small or negative. 
At T > 2 GK, however, the errors for all nuclei with mass number A > 56 are 
below 5% already and thus negligible. Figure U shows target nuclei for a capture 
where errors maximally reach 5-50% at T = 1 and 5-10% ai T — 2 GK, respec- 
tively, when computing the photodisintegration rate from ([28]). Again, this occurs 
for some (but not all, depending on the energy-dependence of the effective cross 
section) a captures with small or negative Q-values. At T = 2 GK even fewer rates 
are affected and the maximal error is below 10%. The shown rates have little or 
no astrophysical impact because those nuclei can only be reached at higher tem- 
peratures. For instance, the 7-process (p-process) significantly photodisintegrates 
nuclei close to stability at T > 2.5 GK and the rp-process als o requ ires such high 
temperatures and probably does not proceed beyond A « lio! ^ ^ 
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Fig. 4. Target nuclei for which an error of 5-50% is introduced in the relation (a, 7)7(7, a) at T = 1 
and T = 2 GK with the standard approximation of the Planck distribution. Maximum errors only 
reach 10% at T = 2. 
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Fig. 5. Target nuclei for which an error of 5-10% is introduced in the relation (n, 7)7(7, n) for 
T > 0.1 and T > 1 GK, respectively, with the standard approximation of the Planck distribution. 



Figure 2] shows target nuclei for neutron capture where errors reach 5-10% at 
any plasma temperature. As expected, captures with low Q-value, either at the 
dripline or in the vicinity of closed shells, exhibit the largest errors but those do 
not exceed 10%. Nuclei so far from stability are expected to be synthesized only at 
higher temperature. For T > 1 GK a smaller number of nuclei shows errors above 
4%. In fact, for most reactions across the chart the error introduced in the reverse 
rate due to the approximation of the Planck distribution is less than 1%. 

Having assured that the approximation of the Planck distribution in the deriva- 



June 10, 2011 0:11 WSPC/INSTRUCTION FILE rauscher 



The Path to Improved Reaction Rates for Astrophysics 19 



tion of does not introduce a considerable error in the obtained rates, it is 
important to realize that further conditions have to be fulfilled to allow the ap- 
plication of the relations discussed above. We started from introducing a thermal 
population 7^^ in p5|) above and consequently all derivations up to here depend on 
the assumption that the excited states in all participating nuclei are occupied ac- 
cording to this population factor. This is a valid assumption for most astrophysical 
plasmas and most nuclei reach thermal equilibrium very rapidly through collisions 
and interactions with photons and other plasma components. However, there are 
some nuclei whic h exhibit long-li ved isomeric states (well-known examples are ^^Al, 
176LU, and i80Taj ^7 | 18 | 19 | 2() | 21 J 

with such spins that they cannot be easily excited 
or de-excited through electromagnetic transitions. At sufficiently high temperature 
they may still get into equilibrium, sometimes through couplings to intermediate 
states, but the relevant transitions have to be carefully studied. This can be achieved 
through an internal reaction network, not connecting different nuclei but rather in- 
cluding the different levels within one nucleus .1^^^ Levels not being in thermal 
equilibrium can be included in regular networks in such a manner as if they were 
a different nucleus. Also in this case, however, the populating and depopulating 
reactions have to be known explicitly. This also applies when ensembles of excited 
states are in equilibrium but the different ensembles within a nucleus are not. Then 
each ensemble can be treated as a separate species in a reaction network and the 
reactions connecting the ensembles have to be included explicitly. 

Although not discussed in further detail here, it is worth mentioning that also 
weak interactions are affected by the thermal population of excited states. For 
example, the /?-decay half-life T1/2 of a nucleus will be changed relative to its ground 
state half-life when the ground state becomes depopulated and excited states with 
different decay half-lives are populated. Thus, the decay "constant" Lp — ln(2)/ri/2 
is actually temperature-dependent 



where = 1/t^ and Tg is the decay lifetime of the excited state. Similar consider- 
ations also apply to other processes, such as electron capture and neutrino-induced 
reactions. 

4.2. Reaction equilibria 

Using the reciprocity relations derived above allows to simplify the full reaction 
networks as defined in ([TJ. Such simplifications are instructive because they en- 
able us to study nucleosynthesis properties which are independent of details in the 
hydrodynamic evolution of the system or even, as shown below, independent of 
individual reaction rates. They usually go along with a restriction to only types of 
reactions in the network which are actually necessary instead of blindly evolving a 
large system of differential equations. Such an approach is not always feasible but 
considerable understanding of nucleosynthesis has been gained in the past through 




(29) 
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such means by circumventing the necessity of computationahy intensive calcula- 
tions. Such simplifications remain important today because it is still impossible 
to couple multi-dimensional hydrodynamic simulations to full reaction networks. 
Furthermore, restriction to the essential often provides a much better insight into 
the physical processes than a brute-force full network calculation. Here we are con- 
cerned with the simplifications because it has to be understood when it is necessary 
to know astrophysical reaction rates and when not. 

Setting the abundance change 1^ = on the left-hand side of ^ implies that the 
sum of all rates destroying the nuclear species i is exactly balanced by all production 
rates and the net change in abundance is therefore Zero, leaving the abundance 
constant. This is called steady flow equilibrium. It is especially useful with reaction 
chains where most rates (perhaps except one) are in steady flow equilibrium. Then 
the slowest reaction sets the timescale of the reaction fiow and all other reactions 
adjust. As long as steady fiow is upheld, no full reaction network has to be solved. 
Rather, the ratios of the steady flow abundances of the involved nuclei are related 
by the ratios of their net destruction rates (or, equivalently, production rates as 
these have to be the same). For illustration, let us assume a chain of reactions 
A^B^C^D^... connecting nuclei through reactions with the same 
projectile, where all net reactions are in steady flow and therefore the same, except 
for the one starting at nucleus A. Then 



and 



(30) 



The slowest rate sets the abundance of B through Ya = —Yb — —rA^B/{Pp\a.NA)- 
The use of a complete set of coupled differential equations is not required anymore 
but the important rates still have to be known. Steady state considerations are 
helpful when investigating hydr ostatic hydrogen burning of stars through the pp- 
chains and the CNO cycles! I I In the past they have also been used for se que nces 
of neutron captures in the s-process on nuclei in between magic numbers.'^ The 
fact that separate steady flows can be assigned to each mass region between closed 
shells has been termed local approximation in s-process studies. 

A slightly different concept is to assume equilibrium between a forward and its 
reverse rate. This is the case when the two rates are equal or, in practice, very close. 
Since stellar rates obey the simple reciprocity relations (P7|) and (1^ . respectively, 
it is trivial to show that 

YAYa _ go9aGQ fmAa\^^ ^-QAa/ikT) ^-32) 



YpYt g^gt \ rngb J 
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for a reaction A + a ^ F + b and 



YAYa _ g^ga Gi 



( 



ruAakT 



) 



3/2 



Yc g^ G§ 



^-QAa/ikT) 



(33) 



for a reaction A + a -H- C + 7. The individual rates do not appear anymore in the 
relation between the abundances. Note that this does not imply that the abundances 
remain constant, they still depend on T which may vary with time as well as on Ya 
and Yf). 

Depending on the plasma densitjQ, above T w 4 — 5 GK all reactions (with 
the exception of the weak interaction) achieve equilibrium. It can be shown that 
the equilibrium abundance of a nucleus A can be calculated from a set of three 
equations 



where A — N+Z is the mass number, Y^, Yp are the abundances of the free neutrons 
and protons, respectively, and m„ the nuclear mass unit. The binding energy of the 
nucleus with neutron number N and proton number Z is denoted by Ba- The 
sums run over all species of nuclei in the plasma, including neutrons and protons. 
Equation ([55]) expresses mass conservation and (155)1 is the charge conservation. 
The unknown abundances Ya, y^, and Yp are obtained with the above equation 
set. Note that reactions mediated by the weak interaction are not included in the 
equilibrium and Yg may be time-dependent. Again, individual rates are not required 
to determine the abundances. 

When all abundances in the network obey the above relations, full nuclear sta- 
tistical equilibrium (NSE) is achieved. In this case, no reaction rates have to be 
known. In realistic cases, more or less extended groups of nuclei are in statisti- 
cal equilibrium and the relative abundances within a group can be described by 
equations similar to p4|) . The different groups are connected by comparatively 
slow reactions not being in equilibrium, which determine the abundance level of 
one group with respect to another group similar to what was shown above for the 
steady state equilibrium. The rates of these slow, connecting reactions have to be 
known explicitly. This is called quasi- statistical equilibrium (QSE). It appears in 
various kinds of high-temperature burning, such as hydrostatic oxygen and silicon 
burning in massive stars and different explosive scenarios. 

^The rate of the triple-a reaction a-\-a-\- a ^-"^^C is very sensitive to the density and will not get 
into equilibrium for Ppi^ < 10^ g/cm'^. 




(34) 




(35) 




(36) 
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A special kind of QSE is the (n,7) — (7,n ) equi librium or waiting point approxi- 
mation, often used in r-process calculations PSEZi This is nothing else than a QSE 
within an isotopic chain, where neutron captures and (7,n) reactions are in equi- 
librium under very neutron-rich conditions (rin > 10^° cm~^) and T k, 1 — 2 GK. 
For the r-process the network is reduced to neutron captures and their inverse re- 
actions, and to decays (with possible subsequent neutron emission). The decays 
are not in equilibrium and determine the timescale with which matter is processed 
from small Z to the heaviest nuclei. When (n,7)— (7,n) equilibrium is achieved, the 
abundances of nuclei within an isotopic chain are connected by 

which connects the abundance of nucleus A (mass number A) with the one of 
nucleus A' (mass number ^-1-1). There is an exponential dependence on the neutron 
separation energy = Q^^ of A'. Also in this type of equilibrium there is no 
dependence on the individual capture or photodisintegration rates. The r-process 
flow to higher elements, however, depends on the /3~-decay rates which connect 
the isotopic chains and are not in equilibrium. They are very slow compared to the 
rates in equilibrium and that is why "waiting points" are established, which are just 
the nuclei (usually only one or two within a chain) with the highest abundances 
according to ((37)) . The r-process cannot proceed until they decay and their decay 
rates have to be known0 

A similar equilibrium, but between proton captures and (7,p), is reached in the 
late phase of the rp-process on the surface of mass accreting neutron stars.l^ There, 
the waiting points are established close to the proton dripline. 

Nucleosynthesis under extreme conditions, such as encountered in some explo- 
sive scenarios, involves exotic nuclei far from stability. According to the above, 
reaction rates are not needed for all of them because reaction equilibria are es- 
tablished at such extreme conditions. Required are nuclear masses (to determine 
Q- values or binding energies) as well as spectroscopic information and nuclear level 
densities (entering the calculation of the partition functions G) . 

But although NSE, QSE, and the waiting point approximation do not contain 
the rates explicitly, they implicitly depend on them because they determine whether 
nuclei are participating in the equilibrium or not. The higher the rates, the lower the 
temperature at which equilibrium is reached. With a time-dependent T evolution, 
this means that the rates determine whether equilibrium is reached earlier with 
increasing T or the freeze-out happens later with decreasing T. 



''The waiting point(s) in two neighboring isotopic chains do not have to be contiguous and therefore 
the notion of an r-process "path" similar to the s-process path is not valid. 
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4.3. Stellar cross sections and experiments 

In principle, stellar cross sections a* as defined in (fTS]) correspond to physically 
measurable quantities, contrary to the purely theoretical effective cross sections 
introduced in ([T9l) . In practice, cross sections tr — cr^=° measured in terrestrial 
laboratories do not include thermal effects. Therefore the rates derived from them 
do not obey reciprocity relations. To derive astrophysical rates appropriate for the 
utilization in reaction networks, laboratory cross sections almost always have to be 
supplemented by theory to account for the additional transitions not included in 
the measurement. 

The energy- and temperature-dependent stellar enhancement factor for the cross 
sections (again obtained via theory) 

fsME,T) = ^^^^^ (38) 

compares the stellar cross section including reactions from thermally populated 
excited states to the cross section obtained with reactions proceeding from the 
ground state of the target nucleus only. Another definition of a stellar enhancement 
factor involves the rates or reactivities, 

/SEFIJJ- ^^^^ Go{T)j;;°aEe-E/ikT)dE ' ^ ' 

In fact, /sEp is the astrophysically interesting quantity because it shows the de- 
viation introduced in the rate at a chosen stellar temperature T when using the 
laboratory cross section instead of the stellar cross section. Then the fraction of the 
ground state contribution to the stellar rate is X — l/(Go/gEp). 

On the other hand, /ggp is supposed to provide information for the experimen- 
talist on how much off the measured a is in comparison to a* at each energy for 
a given stellar temperature. To this end it was quoted in literature occasionally. 
However, a much more useful definition is 

which weights the excited states appropriately for a straightforward comparison. 
This is the same as /spp only if /gfp is independent of E. This is not necessarily 
fulfilled because the energy-dependence of laboratory and effective cross section 
may be different and this would lead to a different evaluation of the reaction rate 
integral. (It turns out that in practice /gpp is often more slowly varying with energy 
than the cross section across the relevant energy window (see Sec. [3]) and that it can 
be approximated by an energy- independent factor in this case.) It is recommended 
to use only either /gfp or /gpp, depending on whether an energy-dependent measure 
is desired or one independent of interaction energy. As has become clear from the 
derivation of the stellar rate in Sec. 14.11 there is no simple relation between Jsef 
and flip. 
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Fig. 6. Stellar enhancement factors /g|;'p for neutron capture on ^*^Os and ^**Os as a function 
of stellar temperature (Tg is in GK). The reciprocals of the normalized partition functions for 
187,188 Qg j^].g j^jgQ shown. 




Fig. 7. Stellar enhancement factors /ggp for ^*^Os(n,a)^*^W as a function of stellar temperature 
(Tg is in GK). 



It should be noted that while the stellar population factors T'^ are normalized 
to Unity and the normalized partition functions Gq cannot become smaller than 
Unity, the three types of stellar enhancement factors defined above can assume any 
positive value, larger or smaller than Unity. 

In which case do we have to expect large deviations 

V = exp (|ln {{avY) - In (M)|) = exp (|ln [fl^f) |) (41) 

from the laboratory value? (The above definition assures I? > 1; without thermal 
effects V ~ 1.) Again, a scrutiny of the effective cross section appearing in ([55]) 
helps to understand the various dependences. 

Let us start with the dependence on stellar temperature. Naively, one would 
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assume that the higher the temperature, the larger the stellar enhancement will 
become. Indeed, this is often the case but there are further intricacies. There are 
two dependences on T appearing in (p9)) . one in the rate in the numerator and 
the other in the normalized partition function in the denominator. The latter is 
monotonically increasing with increasing T. The integration of the product of the 
effective cross section and the MBD may show a different T-dependence, however, 
which even may not prove to be monotonic. This can be understood by combining 
the knowledge of the relevant energy windows from Sec. |3] with definition ()19|) of 
the effective cross section. Larger T shifts the energy window to higher relative 
energy, both in the entrance channel A + a with the relevant Eq becoming larger 
and in the exit channel F + b (or C + 7) with the relevant Eq = Eq + Qau 
becoming larger (see Fig. [2|) . Since the effective cross section sums over transitions 
with relative energies < E^ < Eq and < E^ < Eq , it becomes obvious that the 
higher the relevant energy window, the more transitions are included. This does not 
necessarily result in an increased rate although it often will. It is also conceivable 
that the additionally included transitions at large T have small cross sections and do 
not provide a considerable increase in the rate. Depending on the type of reaction, 
cross sections for transitions already included at small T may also decrease with 
increasing relative energy. In these cases, /glp will be reduced at larger T because 
Go is always increasing. The stellar "enhancement" may even become smaller than 
Unity and not live up to its name anymore. This behavior should not be viewed 
as monotonic, either, because as additional transitions become accessible at even 
larger T the stellar rate again may increase faster than Gq. 

It is to be expected, however, that J^ef decreasing to very small values for 
very large T after having reached a maximum value, i.e., for very high (but not 
necessarily astrophysically important) temperatures I/Gq < /sef ^ ^ ^^^^ always 
be achieved. This is because with increasing relevant energy Eq ' other reaction 
channels become increasingly important, reducing the cross sections of the individ- 
ual transitions. Additionally, more unbound levels will be included in the effective 
cross section. These may lose particles to other reaction channels (e.g., through pre- 
equilibrium emission; see also Sec. 15.3.21) and also not contribute to the effective 
cross section anymore. It was suggested in Ref . fT3l [29l to include only bound states 
in the definition of the effective cross section. Although this may be a good approx- 
imation, it may neglect some transitions which still can contribute to the effective 
cross section at high temperature. Note, however, that unbound states may not be 
in thermal equilibrium although the timescale for reaching equilibrium under stellar 
burning conditions is short compared to the one of a nuclear reaction.'22l A more 
suitable cutoff, if required, would be the energy -Emax appearing as cutoff in the 
calculation of the partition function in ()16|) , although this may include already too 
many transitions with negligible cross sections. Examples for the above considera- 
tions are shown in Figs. El [7] for neutron-induced reactions on Os isotopes. Neutron 
capture on ^^^Os is important to understand the Re-Os cosmochronometer which 
can be used for age determinations in our Galaxy.^SSEQl Modern measurements of 
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Fig. 8. Example for relative effective weights W for levels at excitation energy i?exc and zero 
spin. The weights are independent of stellar temperature but their slopes depend on the maximal 
transition energy as indicated by the labels. 



neutron capture in the astrophysically relevant energy range have reached a preci- 
sion that requ ires th e inclusion of the /g^p correction, even when it is only a few 
tens of percent .1211221 Due to low-lying levels (J'" = 3/2" at 9.75 and 74.3 keV, 5/2" 
at 75 keV, and 7/2" at 100 keV) in isrQs, is higher for ^^'^Os{n,-if^^Os at low 
temperature than for neutron captures on neighboring isotopes. The enhancement 
factors of the reactions ^^^Os(n,7) and ^*^0s(n,7) both rise to a maximum and 
decline from there, soon reaching values below Unity and thus not being "enhance- 
ments" anymore. While the /glp for both capture reactions stay below a factor of 
two (and reach only a few tens of percent at the s-process temperature relevant 
for the cosmochronometer), the /gpp of ^^^Os(n,a)^®^W is also shown here to give 
an example for larger values (see below for the even larger factors encountered for 
photodisintegration reactions). 

Having realized the importance of constraining by the relevant energy win- 
dow at a given temperature and the definition of the range of energies < < Eq 
of relevant excited states, we can arrive at a more general understanding of when 
V will considerably differ from Unity. In the assessment of which excited states in 
the target are contributing to the stellar rate, it is incorrect to directly use the 
Boltzmann population factors from Eq. (|15p ! These are only appropriate when 
each state is bombarded by its own MBD of projectiles. They cannot be used when 
we calculate the rate as usual, integrating over just one MBD with the energy scale 
being relative to the ground state. 

As shown in Eqs. ([23|) -([25 |l the transformation to a single MBD is possible 
because the Boltzmann weight exp{— E'^ / (kT)) in the population factor offsets the 
exponentials in the MBD for each state. This offset is not a complete cancellation 
but results in a transformation of the population P^. The individual weights are 
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transformed to linearly declining functions of the excitation energy 

This is also readily seen in the weighting of the reactions from excited states in 
the definition of the effective cross section by Eq. (|19p . These relative effective 
weights are the ones to be employed when using the standard definition of 
the rate as in instead of the relative weights P^. Contrary to the P^, the 
relative effective weights are not explicitly temperature dependent anymore. 
However, they depend on the energy Eq of the transitions to the ground state which 
are the transitions with the highest possible relative energy. In the application 
to astrophysical reaction rates the range of Eq is given by the relevant energy 
windows as discussed in Sec. [3l It is a relatively narrow range, depending on stellar 
temperature and the type of reaction (and thus introducing an implicit temperature 
dependence). Due to the linearity in excitation energy E^ - contrary to the which 
show an exponential decrease - almost all states up to Eq contribute. Neglecting 
the spin weights g^, one has to consider levels up to E^ « (2/3)i?o to include 90% 
of the levels with non-negligible weight. Depending on the nucleus and considered 
reactions, some levels with large spin values and/or large reaction cross sections 
may require the inclusion of all levels up to very close to E"^ « Eq. The relative 
effective weights as the product of the population factors with their respective 
MBDs are shown in Fig. [H 

Understanding the we are ready to make some general statements on the 
magnitude of 2? in various reactions in a range of astrophysically relevant tem- 
peratures. For example, neutron captures have their relevant energies around the 
maximum of the MBD,Eo « kT. This is between a few keV up to several tens of 
keV for the s-proces&i^ (see also Sec. [3]). Therefore only (exothermic) captures on 
nuclei with low-lying excited states below several tens of keV will exhibit V > 1. For 
light and intermediate mass nuclei, the average level spacings typically are larger 
than « 100 keV and thus V remains close to Unity, with a few notable exceptions. 
The r-process involves neutron captures at temperatures of 1 — 2 GK which 
translates to a location of the MBD peak at « 80 — 160 keV and we will expect 
slightly larger V than for the s-process on average]^ 

Reactions having their relevant energy window determined by charged particle 
widths have considerably larger Eq, in the range of about 0.5 — 13 MeV, depending 
on the stellar temperature and the charges of the interacting particles. Accordingly, 

Since the s-process proceeds along the line of stability, the level structure of the involved nuclei, 
including isomers, is known quite well. This is problematic in the r-process as the experimental 
nuclear structure information far off stability is limited. 
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Fig. 9. Schematic view of the relative effective weights Wa and Wp in the initial and final nucleus, 
respectively, under the assumption of a positive reaction Q- value Q = QAa- In both nuclei a few 
low-lying states are explicitly numbered, above them nuclear level densities are indicated (shaded 
areas). 



T) can already become large at rather low temperature because transitions from 
many excited states have to be considered. Only in light, strongly bound nuclei 
with level spacings of several MeV above the ground state, T) may still remain close 
to Unity. 

The relative effective weights W derived in (|42)) also allow to understand a rule- 
of-thumb which has been known, but never quantified, for quite some time (see, 
e.g., Ref. [13]). The rule states that exothermic reactions (Q = QAa > 0) usually 
have smaller V than their endothermic inverses. Defining the forward reactions by 
the exothermic reaction and the reverse reaction by its endothermic counterpart, 
this means that V""" > V^°™ and frequently 2?'''=^ > pfo™. This is not immediately 
comprehensible upon inspection of the Boltzmann weights as these seem to act 
similarly in the initial and final nuclei. As stated above, however, it is a mistake 
to straightforwardly use the together with a single MBD. The rule becomes 
obvious when using the appropriate W^. The situation is sketched in Fig. [9] where 
the behavior of the weights is shown for two nuclei A and F being the target and 
the final nucleus of a reaction, respectively. For the relevant interaction energy 
set to Eq the maximum relative transition energy is in the target nucleus and 
Eq = Eq+Q in the final nucleus. Provided that Q > the energy range of possible 
transitions < E^, < Eq in the final nucleus is larger than in the target nucleus. 
In consequence, the relative effective weights Wi, in nucleus F decline slower with 
increasing excitation energy and thus a larger range of levels contributes. This 
assumes similar spin structure in the two nuclei, of course, and the rule may not 
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Table 1. Stellar enhan cem ent factors in (7,n) reactions. Values are from a NON- 
SMOKER calculationpH as reported in Ref. |35| 



Nucleus: 


rratc 
JSEF 


Nucleus: 


rrate 
JSEF 


Nucleus: 


rratc 
JSEF 


Nucleus: 


jratc 
JSEP 


186 


400 


190 p^. 


5500 


19Uu: 


1100 


204Hg. 


43 


185 Re: 


1300 


I92pt: 


3300 


196 Hg; 


1700 


204pb. 


160 


l87Re: 


1200 


198pt: 


310 


i98Hg: 


750 







work at small \Q\ when W^^, Wi, are similar and the spin factors g^, gi, dominate. 

The largest reaction Q-values are encountered in capture reactions. For instance, 
neutron captures close to stability exhibit Q- values of the order of 5 — 13 MeV, for 
highly proton-rich nuclei they can reach 20 — 29 MeV. Similar values are found for 
proton captures around stability and on the neutron- rich side of the nuclear chart. 
In the light of the above it is not surprising that endothermic photodisintegration 
reactions exhibit very large /ggp of the order of several hundreds to thousands as 
illustrated in Table [TJ 

An ex ception to the above rule T)'"^^ > T>^°'^'^^ however, has recently been 
discovered.^^SEZl Although the effective weight may be slowly decreasing with 
increasing excitation energy of a level, the corresponding cross section cr'^ may 
decrease much faster, even exponentially. This is because with increasing the 
relative interaction energy in that channel iJ^ = Eq — E^ is reduced. If the cross 
section cr^ is strongly decreasing with decreasing energy - as it is the case in the 
presence of a Coulomb barrier or at high relative angular momentum - transitions 
on excited states will cease to importantly contribute to the effective cross sec- 
tion even when being strongly weighted. This is the reason why charged particle 
reactions show only moderate values of V at low T even though their relevant en- 
ergy window may be high above the reaction threshold. This, of course, acts in the 
entrance channel of a reaction as well as in its exit channel. 

The point, however, is to realize that the barriers may be different in the two 
channels, leading to a different suppression of contributions. For instance, a (n,p) or 
(n,a) reaction has a Coulomb barrier only in the exit channel, {p,a) reactions have 
different barriers in entrance and exit channel. This has the consequence that low- 
energy transitions in the exit channel of an exothermic reaction may be suppressed 
in such a manner as to yield 'D''°^ < Whether this is the case will strongly 

depend on the Q-value because it determines the range of transition energies to be 
efficiently suppressed. The higher the barrier, the larger Q- value is allowed while still 
permitting to suppress most transitions from excited states. The stellar temperature 
also plays a role but its impact is smaller. 

A global search across the chart of nuclides for reactions exhibiting J)""^^ < X>^°™ 
was performed in Refs. I36( 1371 The study focused on identifying cases which are 
interesting for experiments by requiring 'D'^'^^ < 1.5 and the stellar enhancement 
factors for forward and reverse reaction to differ by more than 10%. At stellar 
temperatures T < 4.5 GK about 1200 reactions were found but not all of them are 
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Fig. 10. Reaction Q-values of endothermic (p,n) and {a,ii) reactions which still permit a consid- 
erable Coulomb suppression of the stellar enhancement factors. 

astrophysically interesting.13^ Figure [10] shows the dependence of the Q-values on 
the charge of the target nucleus of two types of reactions selected from the total 
set. The envelope from the maximal \Q\ appearing for each type at each charge 
Z illustrates the action of the increasing Coulomb barrier. It can be clearly seen 
that larger maximal \Q\ is allowed with increasing charge Z. The increase with 
increasing charge is steeper for (a,n) reactions than for (p,n) reactions due to the 
higher Coulomb barrier for a particles. The scatter at a fixed charge number Z is 
due to the range of Q-values found within an isotopic chain. 

Summarizing the above, comparing the position of the relevant energy window 
with the average level spacing in a nucleus already gives an estimate of the mag- 
nitude of the stellar enhancement to be expected. When measuring reaction cross 
sections for astrophysics it is desireable to perform the experiment in the reaction 
direction showing smallest T) in order to stay as close as possible to the actually 
required stellar cross section. This favored direction turns out to be the exothermic 
one in the vast majority of cases, i.e., a reaction with positive Q- value, with com- 
paratively few exceptions (compared to the total number of conceivable reactions) 
as discussed above. It should be remembered that the reverse rate can always be 
calculated using the reciprocity relations (|27|) and (f28| l^ 

"^There is another reason why the knowledge of the exothermic rate is favored which is not con- 
nected to the stellar enhancement. Reaction network calculations implement rates either through 
tables or by using fits of the rates as functions of temperature. In the latter case it is highly advan- 
tageous to be able to perform the fit of the exothermic rate. This is because any small deviation 
of the fit from the actual rate will be strongly enhanced by the factor exp {—Q/{kT)), appearing 
in (|28)l, when Q < 0. 
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Finally, it is worth mentioning that similar considerations for the transforma- 
tion of effective weights and deviations V from the ground state values also apply 
to other types of reactions, not only the ones discussed here. For instance, also for 
electron captures and other reactions mediated by the weak interaction an effective 
cross section, effective weights, and relevant energy windows can be derived, ap- 
plying the appropriate energy distributions. They may become more complicated, 
however, than the ones obtained for the MBD because an explicit dependence on 
the chemical potential is required in the general case. Moreover, neutrinos are not 
in thermal equilibrium with nuclei and thus their temperature will enter as an 
additional parameter. 

4.4. Electronic plasma effects 

The plasma temperature in the astrophysical sites where nuclear reactions occur is 
so high that the nuclei are fully ionized and embedded in a cloud of free electrons. 
This situation affects reactions and decays and has to be considered when preparing 
a rate to be used in astrophysical reaction networks. Similar to the treatment of 
stellar enhancement factors, these corrections also have to be modeled theoretically 
and experimental rates have to be corrected for them. For completeness, some of 
the effects are discussed briefly without going into detail here. 

Decays (and electron captures) are affected not only by the thermal excitation 
of the nucei as shown in (1^ but also by the electrons surrounding the nuclei. A 
nucleus in an atom can be converted by, e.g., capture of an electron from the atomic 
K-shell. In the plasma, electrons are captured from the electron cloud and this may 
alter the half-life of a nucleus considerably. A well-known example for this is the 
decay of ^Be. Its lifetime under central solar conditions (tcc = 140 d) is almost 
double the one under terrestrial, non- ionized conditions (tcc = 77 d). Also the 
decay lifetimes are modified by a change in the electron emission probability. The 
plasma electrons reduce the phasespace available for emission to the continuum and 
thus increase . On the other hand, bound-state decay, i.e. the placement of the 
emitted electron into a low-lying atomic shell, becomes possible even when it were 
forbidden in an atom because of the occupation of available electron shells. Similar 
considerations apply to charged-current reactions with electron neutrinos. These 
effects act in addition to the alteration of the lifetimes through thermal excitation 
of the nucleus and depend on the plasma temperature, density, and composition, 
which may all affect the distribution of the electrons throughout the plasma. 

Nuclear reactions with charged nuclei are affected by electron screening. Elec- 
trons in the vicinity of the nucleus shield part of its charge and thus effectively 
lower its Coulomb barrier. Theoretical predictions of cross sections and rates al- 
ways assume bare nuclei without any electrons and therefore have to be corrected 
for screening. The magnitude of the screening is strongly dependent on the temper- 
ature T, density ppia, and composition Y = J^i Yi of the plasma. At high density, 
electron screening can increase the rate by several orders of magnitude and lead to 
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pycnonuclear burning at much lower temperatures than otherwise needed to ignite 
nuclear burning. Moreover, dynamical effects due to the fast movement of electrons 
may also play a role, although there is discussioi i on whether this is important when 
all reactants are in thermodynamic equilibrium It is apparent that the theoreti- 
cal treatme nt of scree ning is complicated and we are far from complete microscopic 
descriptions P^l^nSIlTJig proper inclusion of screening effects remains a challenging 
problem in plasma physics. 

Under most conditions, the screened reactivity can be decomposed into the 
regular stellar reactivity and a screening factoJ" 



screened 



:C(r,ppia,r)M* . (43) 



This implies that the Coulomb potential Vc seen by the reaction partners can be 
described by the bare Coulomb potential and an effective screening potential U 

Vc{r, T, ppia, Y) = V^-'%r) + t/(r, T, p^u, Y) = + C/(r, T, p^u, Y) . (44) 

Then the screening factor acquires the form C = exp(— {/ /{kT)). The challenge is 
in the determination of U depending on the plasma conditions. 

An often used static approximation, being appro priate for early burning stages 
of stars, is weak screening in the Debye-Hiickel model SSlll Weak screening assumes 
that the average Coulomb energy of each nucleus is much smaller than the thermal 
energy, i.e. ZeVc <SC kT. Then a nucleus will be surrounded by a polarized sphere 
of charges, with a radius 



with 



/ kT 
y Aiie^ppuNAC 



C = 5^(Zf + Z0V, , (46) 



where the sum runs over all charged plasma components. The screening factor C is 
transformed to Co = 1 - UT)/{kT) with C/d = -e^Z^Za/iiD- 

For strong screening in high density plasmas it is more appropriate to use the 
ion-sphere model instead of the Debye-Hiickel approximation.^^ The ion-sphere 
model is equivalent to the Wigner-Seitz model used in condensed matter theory. 

Another type of screening is observed in nuclear experiments in the laboratory. 
There, nuclei are present in atoms, molecules, or metals, each with specific electron 
charge distributions around the nucleus. Although completely different from plasma 
screening, this type of screening has to be understood because it is especially im- 
portant at the low interaction energies of astrophysical relevance. The measured 
reaction cross sections have to be corrected to obtain the bare cross section which 
can be compared to theory or used to determine the rate. Atomic screening can be 
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treated in the adiabatic approximation, leading to 



■ exp —— (47) 



a{E) E + Ua \ E 

with the Sommerfeld parameter 77 from (|lip . The screening potential Uc in this 
approximation is given by the difference in the electron binding energy of the target 
atom and the atom made from target atom plus projectile Ue = B^~^°' — . In light 
systems, the velocity of the atomic electrons is comparable to the relative motion 
between the nuclei. Therefore a dynamical model is more appropriate!^ However, 
the adiabatic approximation provides an upper limit on the expected screening 
effect on the cross section. 

There seem to be discrepancies between theory and laboratory determinations 
of Uc, the latter often yielding much larger values of Ug. Some of them have been 
resolved through improved stopping powers used in the determination of the exper- 
imental cross sections, '^^ ^'^ while others r emain puzzling , especially regarding 
cross sections of nuclei implanted in metals.^^ ^ l ^^ l ^^ l ^^ l ^^ l Thus, the laboratory 
screening seems to be less understood than the stellar screening. 



5. Reaction mechanisms 
5.1. General considerations 

Having discussed the special requirements of astrophysics for the determination of 
the stellar rates in the preceding sections, we turn to the question of how to ob- 
tain the cross sections a^'^ required in Eqs. ([8]), (fT5|) and p9)) at the energies of 
astrophysical relevance. It has become apparent that astrophysical rates include 
more transitions than usually obtained in straightforward laboratory reaction cross 
section measurements. Cleverly designed experiments may study some of them but 
especially at larger stellar temperatures, as they are typical for explosive burning, 
theory will be indispensable for providing the appropriate stellar reaction rates. 
Even at stability, the small cross sections of charged-particle reactions at astro- 
physically relevant energies (see Sec. [S]) pose a considerable challenge for current 
measurements and future experiments have to employ novel techniques or new fa- 
cilities to address this problem. Moreover, hot astrophysical environments produce 
highly unstable nuclei which cannot be studied in the laboratory, yet. Reaction 
cross sections of nuclei far from stability at astrophysical energies probably will 
never be experimentally determined. Therefore reaction networks for explosive nu- 
cleosynthesis have to include the majority of their reaction rates from theoretical 
predictions although experiments may help to determine nuclear properties and 
some of the transitions required for the calculation of effective cross sections. 

Although reaction theory dates back as far as the 1950s, the special requirements 
of astrophysics and the need for cross section predictions of nuclei far from stability 
provide an interesting and stimulating environment for the application and further 
developments of different approaches. The challenges are manyfold. On one hand. 
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astrophysical energies are very low and, as we will see below, different reaction 
mechanisms may contribute or interfere. On the other hand, even if the reaction 
mechanism is unique and well understood, nuclear properties entering the reaction 
model have to be predicted for nuclei far off stability. This proves challenging even 
for modern nuclear structure calculations. Although a fully microscopic treatment 
is preferrable, good parameterizations and averaged quantities are still necessary 
in many cases due to the sheer number of reactions and involved nuclei, especially 
for intermediate and heavy nuclei consisting of more than 30 — 40 nucleons and 
thus not allowing the application of few-body models. Finally, the interpretation of 
experiments has to be supported by theory. This latter case may involve different 
methods than the one dealing with the prediction of astrophysical rates because 
experiments may be conducted at higher e nerg ies and theory is needed to extract 
the information to be included in the rates.!™ For example, the properties of ex- 
cited states and their spectroscopic factors can be studied by (d,p) reactions at 
comparatively high energy which are not directly relevant in astrophysics. In the 
following I focus on theory for the prediction of astrophysical react ion r ates. 
Theoretical models can be roughly classified in three categories f^^^ 

(1) Models involving adjustable parameters, such as the R-mdXYY^^ or the K- 
matrijI^Sl methods; parameters are fitted to the available experimental data and 
the cross sections are extrapolated down to astrophysical energies. These fitting 
procedures, of course, require the knowledge of data, which are sometimes too 
scarce for a reliable extrapolation. 

(2) "Ab initio" models, where the cross sections are determined from the wave 
functions of the syste m. The potential model^, the Distorted Wave Born Ap- 
proximation and microscopic 

niodel J58|59l60l 

are, in principle, 

independent of experimental data. More realistically, these models depend on 
some physical parameters, such as a nucleus-nucleus or a nucleon-nucleon in- 
teraction which can be reasonably determined from experiment only. The mi- 
croscopic Generator Coordinate Method (GCM) provides a "basic" description 
of a nucleonic system, since the whole information is obtained from a nucleon- 
nucleon interaction. Since this interaction is nearly the same for all light nuclei, 
the predictive power of the GCM is high for such nuclei. 

(3) The above models can be used for low level-density nuclei only. This condition 
is fulfilled in most of the reactions involving light nuclei {A < 20). However 
when the level density near threshold is large (i.e. more than a few levels per 
MeV), statistical models, using averaged optical transmission coefhcients, are 
more suitable (see Sec. I5.3.2t . 

The nuclear level density (NLD) at the compound nucleus excitation energy cor- 
responding to the astrophysical energy window determines which reaction mech- 
anism is applicable and which model to choose. The compound formation energy 
^-form = Eq + Escp IS givcu by the astrophysical energy Eq relative to the ground 
state of the target nucleus and the separation energy of the projectile -Escp- Around 
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stability, Egcp usually is high and dominates -Eform- Statistical models will then be 
applicable for intermediate and heavy nuclei with sufficiently high NLD. Even at 
stability, however, the NLD may not be high enough at nuclei with closed shells 
(see Sec. I5.3.2|) . Approaching the driplines, the neutron- or proton-separation ener- 
gies strongly decrease, resulting in low -Bform for neutron- or proton-induced reac- 
tions, respectively. This leads to low NLD at £'form even for intermediate and hea vy 
nuclei.aisolated' resonances but also direct reactions will become important.'^ 

In the following, reactions at intermediate (Sec. 15. 2|) . high (Sec. 15. 3p . and low 
(Sec. 15. 5p compound NLD are discussed separately although there may be con- 
tributions from several reaction mechanisms simultaneously, especially in systems 
with low and intermediate NLD. The discussion will focus on models more or less 
applicable for large-scale predictions across the nuclear chart. Further models are 
presented in, e.g., Refs. [71 [53 l [59 l l60l 



5.2. Resonant reactions 

Resonances in reaction cross sections are important for the majority of nuclei. 
Depending on the number of nucleons in the target nucleus resonances appear in 
the reaction cross section at lower or higher energy and their average spacing also 
depends on the structure of the nucleus. Astrophysical energy windows cover regions 
of widely spaced, isolated resonances to regions of a large number of overlapping, 
unresolved resonances. Accordingly, different approaches have to be combined. The 
latter region is more important for reactions between charged reactants because the 
Coulomb barrier shifts the relevant energy window to higher energy compared to 
reactions where neutrons determining the location of the window (see Sec. [3]). 

Isolated resonances in the low and intermediate NLD regimes can be treated in 
the i?-matrix^ or the i^-matriJ^SI approaches or by applying simple single-level 
or multi-level Breit-Wigner formulae. In all these methods, the resonance proper- 
ties (resonance energy, spin, partial and total widths) have to be known. Often, 
an inverse approach is used and the resonance properties are derived from experi- 
mental data by, e.g., i?- matrix fits. Where this is impossible, nuclear theory has to 
be invoked to predict the required quantities. This remains problematic, however, 
because the reaction cross sections are very sensitive to the resonance properties. 
The astrophysical rates are also sensitive but since their calculation involves an in- 
tegration over an energy range, only strong resonances truly contribute and others 
may be averaged out. Nevertheless, there are large uncertainties in reaction rates 
off stability due to the unknown resonance contributions (see also Sec. l5.4|) . Cluster 
models (see Sec. 15. 3. H and Ref. l53|) have been successful in describing resonant cross 
sections in light nuclei but cannot be easily applied to nuclei at intermediate and 
heavy mass. 

Although resonances with the same spin J interfere and single resonances may 
also show interference with a direct reaction (Sec. l5.5|) . the single- level Breit-Wigner 
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formula (BWF) is often usedl^UJJ 



It is quoted here despite its restrictions because it allows to demonstrate some 
important principles important for the relation between stellar and laboratory rates. 
Equation (|48p gives the BWF for n non-interfering resonances. The wave number is 
denoted by . The total width r*°* of a resonant state i in the compound nucleus is 
the sum over the widths of the individual decay channels r*°* = r^"+rf'' + . . . , also 
including transitions to other reaction channels beyond the exit channel F + b. The 
widths of the individual decay channels are summed over transitions to all possible 
final states in the channel. Thus, Ff^ = Tlu^'i- Figure [2] shows the energy scheme 
and the contributing transitions in each channel. When the resonance energy E'^"^ 
is known, the widths F and F can be calculated from the transmission coefficients 
obtained by solution of a Schrodinger equation in the optical model (see Sec. 15.3. ip 
and a spectroscopic factor (see Sec. l5.5p . Both resonance energies and spectroscopic 
factors should, in principle, be predictable in the shell model (see, e.g., Refs. IMl 
l63l [64|) or other microscopic theories but this is currently not feasible for all nuclei 
across the nuclear chart. Different approaches yield results which differ more than 
it is tolerable in the calculations of astrophysical reaction rates.^^ 

According to ([TSl) . we take ^ = for the usual laboratory cross section. It is 
interesting to note that also resonances located below the reaction threshold may 
contribute due to their finite width reaching above the threshold. These are called 
sub-threshold resonanceJ^^ In reactions with a large, positive Q- value the energy 
dependence of the partial width in the exit channel Ff can be neglected. This is 
not true when the Q-value is small or negative. 

As has become obvious in the discussion of the stellar cross section in Sec. 14.11 
for the astrophysical reaction rate, the effective cross section has to be employed in 
the integration for the reaction rate and thus a weighted sum over excited target 
states fjL has to be performed and we obtain 



1/2 ^ ^ A 

) {kT)-^/^ j\y2y^,<w{E,)\E^^-^dE^ , (49) 



with the effective weights W taken from p2)) and E^^ — — E^. This can be sim- 
ply achieved by replacing F^ in (^51) by F^" (summing over all possible transitions 
in the entrance channel) and dividing the resulting cross section by the normal- 
ized partition function of the target nucleus Gq. This can easily be shown when 
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combining definitions (jl9p and (j48l) for a single resonance with spin J, 

^ go Eo gg^^^^ kl {E^ - £;res)2 + (ftot/2)2 

.goga^^o ^ 2mAa (£;o - E'^""")^ + (ft°t /2)2 
^ 7r ffj(l + (5Aa) 



This results in 



/BW 



/ 8 



/BW.i 
1/2 



/ Stt 



CO 



1/2 



l+<5Aa 



r^(i?„^)rf^(ii;, 



■(^A_^re.)2 + (rtot(^A)/2)2 



E^e'-^ dE^ 



(51) 



for the stellar rate. 

Some simplifications can be made depending on the resonance widths. For sim- 
plicity, the derivations are given for a single resonance with spin J, for n > 1 
the contributions can be added. A frequently used simplification is the one for 
narrow resonances, assuming that the widths F in the numerator in the integral 
and the Boltzma nn factor exp{EQ/{kT)) do not change across the width of the 
resonance .^2lll29l Then their values can be taken at the resonance energy and the 
integration can be performed analytically (see also Eq. ^59\i ). yielding 



,y ( 27r \^'\ 1 {l + 5Aa)gjr^'^{Er)r^''{Er) 

^ \mAakTj ''^ gago ft°t(i?r) 

(52) 

where the resonance energy Eq"^^ is given relative to the ground state of the target 
nucleus. How does this compare to the ground state rate usually measured in the 



June 10, 2011 0:11 WSPC/INSTRUCTION FILE rauscher 



38 T. Rauscher 



laboratory? The stellar rate can be recast as 



{av): 




with 

\m,AakT J gago T*-°^ {E^'''') 

The second line in ([55)) was obtained by neglecting the energy dependence of T^^ 
and r'°*. This is a vahd assumption provided the reaction has a sizeable, positive Q- 
value. As can be seen from (j53p . the contributions from excited states vanish quickly 
because often F is strongly energy dependent and vanishes fast with decreasing 
energy (remember that E^ — Eq — E^). This is certainly true for reactions between 
charged particles and low resonance energies. Nevertheless, the resonant transitions 
from excited states may dominate a resonant stellar rate.^ 

For broader resonances the above approximation cannot be used and the rate 
has to be determined by numerical integration of Eq. (|5]). The wings of broad 
resonances can also contribute significantly to the rate even when the resonance 
energy is outside the relevant energy window for the rate. Sometimes the values 
for F", F^^, and r*°* are known experimentally at the resonance energy. Then a 
frequently used approach in experimental nuclear physics is to assume that F^'' and 
F'°* are approximated by energy-independent values and the energy-dependence 
of F*^ is only due to a barrier penetration factor derived from an optical model. 
Even if the energy dependence of F^^, r*°* is accounted for explicitly, this type 
of extrapolation does not include the stellar enhancement and is only valid for 
laboratory cross sections. The stellar rate must be calculated from a weighted sum 
of resonant contributions, as shown in ()49|) . both for the value at the resonance 
energy and in the extrapolation. It follows from ()51|) that the only difference in the 
energy dependence, however, stems from the width in the entrance channel where 
F° has to be replaced by F'^". The additional transitions to excited states in the 
target nucleus can be measured in principle. If they are not available, F"^" at the 
resonance energy has to be predicted from theory. Also the extrapolation is more 
involved because F"^" will have a different energy dependence than F°. The same 
methods can be used as in the extrapolation of F*^ but they have to be applied to 
all contributing transitions separately. 

The discussion of electron screening in the stellar plasma in Sec. 14.41 applied 
to nonresonant rates. It can be shown that the same screening corrections can be 
applied for resonant rates when f"^° ^ f^''.'^' A more complicated form arises for 
pAa ^ fFb^ Refs.EllandEHlfor further details. 
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5.3. Resonant reactions at large compound level density 
5.3.1. Optical model 

Microscopic models are based on basic principles of quantum mechanics, such as 
the treatment of ah nucleons, with exact antisymmetrization of the wave functions. 
The hamiltonian of an A-nucleon system is 

A A 

^^E"^ + E ^^^^ (55) 

4=1 i<j = i 

where % is the kinetic energy and Vij a nucleon-nucleon interaction.l^SEl] ^pj^g 
Schrodinger equation associated with this hamihonian can not be solved exactly 
when A > 3. The Quantum Monte Carlo method represents a significant break- 
through in this direction, but is currently limited to A — 10.^^ In addition its 
application to continuum states is not feasible for the moment (it has been applied 
to the d(a,7)^Li reaction but the a+d relative motion is described by a nucleus- 
nucleus potential).!^ 

In cluster models, it is assumed that the nucleons are grouped in clusters and 
internal wave functions describing the relative cluster motions are generated.'^ The 
main advantage of cluster models with respect to other microscopic theories is 
its ability to deal with reactions, as well as with nuclear spectroscopy. Over the 
past years, much work has been d evoted to the improvement of the intern al w ave 
functions: multicluster descriptioni^lEl], large-basis shell model extensionJ^, or 
monopolar distortiorPSl. The main limitation arises from the number of channels 
included in the wave function, which reduces the validity of the model at low ener- 
gies. Also large NLDs require many channels in the wave functions. Therefore the 
application of cluster models is limited to light nuclei. 

Due to the complexity of the nucleon-nucleon (NN) interaction, one often resorts 
to working with effective interactions instead of solving microscopic models based 
on NN potentials. Widely used in calculating different reaction mechanisms is the 
optical model. In that model, the complicated many-body problem posed by the 
interaction of two nuclei is replaced by the much simpler problem of two par ticles 
interacting through an effective potential, the so-called optical potential.''''^ 76|77| 
Such an approach is usually feasible only with few contributing channels. Always 
included is elastic scattering. That is why optical potentials can be derived from 
elastic scattering data. 

The time-independent radial Schrodinger equation is numerically solved with 
an optical potential which provides a mean interaction potential, averaging over 
individual NN interactions. In consequence, single-particle resonances cannot be 
described in such a model. However, resonances stemming from potential scatter- 
ing can still be found. Elementary scattering theory yields expressions for the elastic 
cross section and the reaction cross section. The latter includes all reactions and 
inelastic processes which cause loss of flux from the elastic channel. With the di- 
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agonal element 5'"" of the S-matrix (sometimes also called scattering matrix or 
collision matrix), the reaction cross section for spinless particles is then given bjE3 

<" = ^E(2^+I)(l-I^ri') • (56) 

This can be generalized to other outgoing channels /?, not just the elastic one. The 
elements of the S-matrix are complex, in general, and related to the scattering 
amplitude / of the outgoing wave function 

r'^^Y.^'2^+i){sf-m{cose) , (57) 

which, in turn, is nothing else than the transition amplitude t"'^ — -~{2'Kh)/mpf°'^ 
connecting entrance channel a and exit channel /3, with m being the reduced mass 
in the entrance channel and Pi a Legendre polynomial. The imaginary part of the 
optical potential gives rise to an absorption term in the solution of the Schrodinger 
equation, thus removing flux from the considered channels. Therefore, the matrix 
element S*"" is also related to the transmission coefficient 

Tt=(l- |e2'**f ) = 1 - l^ri' (58) 

which describes the absorption of the projectile by the nucleus. Important for prac- 
tical application is that the phase shifts 5i can be derived from elastic data. 

The optical model is well suited for describing transitions between states of 
intermediate and heavy nuclei. It has been and is still used also to treat reactions 
with light nuclei although other methods exist for these. 



5.3.2. Statistical model 

The optical model can be used to compute the widths (F, F, 1"*°*) appearing in 
the BWF for resonances, see (|48| above. As mentioned before, the relevant energy 
windows for astrophysics also include compound nucleus excitation energies with 
such high NLD that individual resonances cannot be separated because the average 
resonance width (F) becomes larger than the average level spacing D = 1/p. In 
fact, this is the case for the majority of reactions included in astrophysical reaction 
networks. Instead of explicitly dealing with a large number of unknown resonances, 
one moves to averaged resonance properties. 

Starting with the BWF given in , the sum of individual resonances can be 
replaced by an average over an energy interval A.E using the mathematical relatiorP 

' ' dE,,, 



(E^. - El^p^ + (Ft°V2)2 / AEJ (E^ _ i?-)2 + (rtot/2) 

. (59) 

AE r'°* 
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Here, the angle brackets denote the average as defined by the above equation. Note 
that the approximation for narrow resonances, as also used in (j52[) and (j54[) . was 
applied to arrive at the last line of (|59|) . This is obviously allowed because of the 
assumption of a large number of narrowly spaced resonances. With this we rewrite 
the sum over resonances in the BWF as 

\V (E„-E.,„f + (rr/2)V \ rg / 



Jtt I 



(60) 



The number of resonances Arinuc within an energy interval Ai? was replaced by 
the NLD p in the last line. The averaged widths, the NLD, and the W are energy- 
dependent, of course. The width fluctuation coefficients W account for the different 
averaging in the last line 



\ ^r.iE) I {r;,(i5)) (f«(B)) 

In terms of physics, they describe non-statistical correlations between the widths 
in the channels A + aand F + b. In practice, they differ from Unity only close to 
channel openings .1^^!^ 

Making use of the relation between transmission coefficients obtained from the 
solution of the Schrodinger equation with an optical potential and the averaged 
widths = 27r/9 (V^) , . . . the cross section for the statistical model of compound 
reactions can be written as 

E ..-^I^— . (62) 

The summation in the denominator runs over all channels c leading to the same 
compound nucleus, not only (but including) A + a and F + b. Thus, this sum 
is equivalent to 7}°*. Also, the sums over channel spins j and partial waves £ are 
explicitly written to emphasize that the transmission coefficients must include these 
quantum numbers. Each transmission coefficient includes transitions from states at 
the compound energy -Eform = E + Escp (i?sep being the separation energy of the 
projectile in the compound nucleus). While T^^ only includes those to the state fi in 
the target nucleus, the T include all transitions allowed by energetics and quantum 
selection rules. 

Comparing (|62|) with (|48l) and (|54|) it is readily seen that the statistical model 
cross section is an averaged Breit-Wigner cross section for narrow resonances, when 
W = 1. Completely equivalently to the Breit-Wigner case, it can be shown that for 



June 10, 2011 0:11 WSPC/INSTRUCTION FILE rauscher 



42 T. Rauscher 



the calculation of the stellar rate it is sufficient to replace T'^ by T"^" in (I62 p and 
divide the resulting integral by the normalized partition function Gq , giving 

1 l+5Aa 

X 




E^e'T^ dE^ . (63) 



The total transmission coefficients in each channel z =Aa, Fb, ...include a 
sum over final states ^ in that channel. Similar to the treatment of the partition 
functions in (|16p . the sum over discrete states can be extended by an integration 
over a level density above the energy i?^i.,3t of the last discrete state included, 

f ^(£;, J, ^) - <j ^ r«(s, J, TT, j«, 

J2T''{E,J,TT,E'-,J\Tr')p'{E',r-,TT')dE' . (64) 

The integration is over the NLD in the channel z, i.e. in the target nucleus A 
for channel Aa, in the final nucleus F for channel Fb, and so on. The transmission 

is the same as T^, only that it is a transition to an artifical state with given 
{E^, J^,TT^) and weighted by the NLD p^{E^, J^, tt^). The relative transition energy 
in channel z is E^ = E—S^^p—E^ = Efomi^E^, where S^^p is the channel separation 
energy. The reader is advised to consult Fig. [2] to get an overview of the included 
transitions and their relative energies. 

Particle transmission coefficients have to obey spin selection rules and thus 

T«(i?,J,7r,i?«,J«,7r«)= E ^i^^^D ■ (65) 

M3-s\j = \Ji-J\ 

Here the angular momentum £ and the channel spin j = J + are connected by 
j = £ + s including the particle spin s. Each Tji can be directly obtained from the 
solution of the (time-independent, radial) Schrodinger equation at the energy i?| 
with an appropriate optical potential. 

The calculation of radiative transmission coefficients proceeds equivalently to 
([65| but electromagnetic selection rules (see, e.g.. Appendix B of Ref. .3) have to 
be obeyed. The parities tt, tt^ and the angular momentum £ select the type of 
allowed electromagnetic transition (El, E2, Ml, M2, etc.) and accordingly the ap- 
propriate description has to be invoked to calculate the transition strength Tji. To 
phenomenologically account for pre-equilibrium particle emission at higher com- 
pound excitation energy (see Sec. I5.3.3|) . sometimes the integration in (|64|) is only 
carried out to a cut-off energy min(i?cut, ^form), with an appropriately chosen i?cut 
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(e.g., the energy at which the 7 transmission exceeds a certain fraction of the total 
transmission) for the 7 transmission coefficient appearing in the numerator of (1631) . 
The total transmission coefficient Tj°' in the denominator, however, always has to 
include the full integration up to -Eform- 

The statistical model of compound reactions was initially developed by Bohr, 
who conceived the independence hypothesis.'" It states that the projectile forms 
a compound system with the target, shares its energy among all of the nucleons, 
and finally the compound nucleus decays by emitting photons or particles indepen- 
dently of the formation process. This implicitly requires long reaction timescales 
as the compound nucleus has to live long enough to establish complete statistical 
equilibrium among the nucleons. Compared to the direct mechanism (Sec. 15. 5|) the 
timescale is about 5 — 6 orders of magnitude longer and includes many degrees 
of freedom. In the independence hypothesis, the (laboratory) cross section can be 
factorized into two terms 

Aa—^Fb Aa i _ Aa \ I _ Aa (aa\ 

CTfjF — Crfoim''doc — CTform T — ^form;f~7 ' \^^) 

\y tot/ J 

the formation cross section of'^^^ and a branching ratio describing the probability 
for decay to the observed channel Fh. An early implementation of this was the 
Weisskopf-Ewing theory.l^ Since then, the Hauser-Feshbach approach has been 
widely used, which also incorporates conservation of angular mo men tum partially 
lifting the independence assumption but thus being more realistic. ISni Equation ^ 
is the cross section from the full Hauser-Feshbach formalism. Nevertheless, although 
too simplified, Eq. (1661) is sometimes useful when estimating the relative feeding of 
different reaction channels. 

Although it might seem tempting to conclude that the cross section of a reac- 
tion proceeding through the compound mechanism should be smooth because it is 
formed from the superposition of amplitudes from a very large number of states 
with random phases, this is a wrong assumption.!^ It was first shown in Ref. |81] 
that the cross sections can continue to show large fluctuations. The usual Hauser- 
Feshbach equations do not account for these fluctuations. Therefore, a meaningful 
comparison to experimental data is only possible after averaging the data over a 
sufficiently wide energy range, comparable to the average resonance widths. When 
using the statistical model to compute astrophysical reaction rates (or when de- 
riving rates experimentally directly) this is taken care of automatically. However, 
when using beams with a very narrow energy spre ad it s hould be noted that the 
results cannot be directly compared to calculations .^^^^^ 

It is worthwhile to point out that the reaction rate is rather "forgiving" to 
deviations around a "true" cross section value, provided the deviations go both 
ways and can cancel within the integration in ([8]). Therefore, the statistical model 
approach may even be applicable in the presence of small but isolated resonances 
as long as their average contribution is correctly accounted for. 

This is closely connected to the question of the applicability of the statistical 
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dows for (11,7) reactions have not been changed!^ For comparison, Fig. [TT] shows 
the minimal temperatures required for (p,7) (although on a coarser temperature 
grid) for a region of the nuclear chart when applying the correct energy windows. It 
has to be kept in mind, though, that reactions with different exit channels do not 
exhibit the same relevant energy window, as explained in Sec. [31 The general pic- 
ture arising is the same as described in Ref. IH The statistical model can be applied 
to the majority of neutron capture reactions, with exceptions close to magic neu- 
tron numbers or low neutron separation energy, both leading to a low NLD at the 
compound formation energy -Eform- Since the relevant energy windows for charged- 
particle capture are shifted to higher energy with respect to the ones for (n,7), 
the applicability is even broader. This may also apply to reactions with a charged 
particle in entrance or exit channel but no general statement can be made because 
it depends on which width determines the location of the energy window. Finally, 
endothermic reactions always require higher temperature to have appreciable rates 
but that does not necessarily mean that the compound nucleus is formed at high 
excitation energy. Therefore, they may require even higher temperature until the 
statistical model becomes applicable. 

5.3.3. Modifications of the standard statistical model 

In the preceding section, the applicability of the statistical model (Hauser-Feshbach 
model, HFM) due to the required average NLD has been discussed. Modern reac- 
tion theory knows a multitude of models, each suited to a particular reaction type 
and mechanism. These all imply certain approximations. Of course. Nature is con- 
tinuous, many types of reactions occur simultaneously and we have to choose an 
approximation suited to describe the dominant effects. Astrophysical energy win- 
dows prefer low projectile energies but still may include transitional regions between 
several reaction mechanisms. 

Discussed below are two types of modifications to the rate calculations: 1) Ac- 
counting for additional reaction mechanisms, and 2) modifying the HFM itself to 
provide a smooth transition to low NLD regimes. The inclusion of isospin con- 
servation is a further modification which is discussed in a separate subsection of 
Sec. [5X21 

Direct reactions are known to be relevant at high projectile energies but were 
also found to significantly contribute at low energies in nuclei with low NLD. They 
are further discussed in Sec. 15.51 

Semi-direct^^ and multistep reactioni^ also occur on a faster timescale than 
equilibrated compound reactions. They become important at several to several tens 
of MeV projectile energy and are thus outside the astrophysically relevant energy 
range. The capture of a-particles on intermediate and heavy nuclei may barely reach 
such energies but only at high plasma temperature. It was shown that semi-dir ect 
capture is negligible even for very neutron-rich nuclei at r-process conditions.'^ 
When assessing the importance of additional mechanisms it is essential to consider 
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the change they bring about for the astrophysical rates. Even if the ground state 
transitions (as usuahy explored in theoretical and experimental nuclear physics) 
may barely reach the required high energy for additional processes, such as semi- 
direct reactions, this may not affect the reaction rates much. This is because they 
also include transitions from excited states according to ((25|) and ((42|) . which pro- 
ceed at lower relative energy and thus remain unaffected by modifications of the 
cross sections at higher energy. 

The HEM assumes that the compound nucleus has sufficient time to distribute 
the energy gained through the interaction with the projectile among all nucleons of 
the compound system before it decays. At a high compound formation energy E'form 
of several tens of MeV transitions occurring before this energy-equilibration lead to 
multistep reactions and pre-equilibrium emission of particles and photons.^^ The 
required compound excitation energy is too high to be of astrophysical relevance. 

The second type of modifications introduced above is discussed in the following, 
focussing on two ideas. The first is that subsequent emission of several particles may 
also occur when the final nucleus is produced at an energy above a further parti- 
cle emission threshold. Astrophysically this can be relevant in nuclei close to the 
driplines with low separation energies. It can be treated approximately by applying 
an iterative application of (|62|) , (|64| , and (|65l) to transitions to states in each system 
formed in each emission process. This has been used, for example, for calculating 
neutrino-induced particle emission which is relevant to the construction of neutrino 
detectors. ^'^ In such calculations, the formation transmission coefhcients Tj in 
([S^ are replaced by neutrino transmission coefficients '^'"^<=Tj describing the pop- 
ulation of compound states by neutrino interactions (using, e.g., the random phase 
approximation). Neutrino reactions select high excitation energies and thus multi- 
particle emission will also be important at stability. The same approach can be 
used to determine /3-delayed particle emission when /3-decay produces a daughter 
nucleus in an excited state above a particle-separation energy. It can als o be u sed 
to study /3-delayed fission. Both processes are important in the r-process.l^SF33l 

The second idea revolves around the fact that regular HEM calculations as- 
sume a compound formation probability independent of the compound NLD at the 
compound formation energy. Therefore the sum in (jl4p runs over all J'^ pairs (a 
high-spin cutoff is introduced in practical application of the model because spin 
values far removed from the spins appearing in the initial and final nuclei do not 
contribute significantly to the transmission coefficients). The availa bility of com- 
pound states and doorway states defines the applicability of the HEM.^I^ Relying 
on an average over resonances, the HEM is not applicable with a low NLD at com- 
pound formation. Furthermore, not all spins and parities will be available with equal 
probability at each inform, especially at low NLD. On average the HEM will then 
overpredict the resonant cross section (unless single resonances dominate) because 
it will overestimate the compound formation probability. This can be treated by 
introducing a modification of the formation cross section which includes the com- 
pound NLD dependence. The summands of (fT4|) will then be weighted according 
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to the available number of states with the given J'^. (Formally this is the same as 
assuming dependent potentials for particle channels.) In the most general case 
this will require a dependence on a J'^-dependent compound nucleus NLD at inform- 
Note that the standard HFM only includes the NLD of the compound nucleus in 
the photon transmission coefficients (see Sec. I5.4.2P to determine the endpoints of 
the 7-transitions. 

The assumption that all spins and parities are available can be lifted in sev- 
eral steps. A parity-dependent, global NLD was used in the calculation of 
the transmission coe fficients f or nuclei without experimentally determined excited 
states (see Sec. 15.4.2"]) j ^^ l ^^ l ^Q ^Shortlv thereafter, the parity-dependence of the com- 
pound formation was implemented in a modified HFM in version 4.0w of the code 
NON-SMOKERWEB ^ggg j^igQ [5X3]) .1^ A discussion of the imphcations of a 
parity-dependent compound formation for astrophysical neutron capture is given 
in Ref. [92l Although the parities are not equidistributed up to sizeable excitation 
energies, the impact on stellar rates remains small (in comparison to other un- 
certainties) because the effective cross section ([T9|) also includes transitions from 
excited states which washes out the selectivity on parity. This gives rise to factors 
of two modifications close to the neutron dripline (but see Sec. 14. 2|) . 

Additionally, NON-SMOKERWEB offered the option of weighting the HFM 
cross section by a function depending on the total NLD since version 4.0w.l2iJ An 
improved implementation with J'^ dependent weighting of the summands, thus im- 
plicitly accounting for a low NLD at the compound formation energy, is introduced 
in the SMARAGDcodeand will be used for a future update of large-scale reaction 
rate predictionsl ^^ l ^^ l ^^ l ^^ l Preliminary results with this modification are shown in 
Fig. [331 hi Sec. [5X51 

Obviously, these modifications of the HFM depend on the NLD treatment to 
obtain the spin- and parity distribution. See the paragraphs on the NLD in Sec. 15.4.^ 
for further details. 

5.4. Sensitivities of HFM rates to nuclear properties and other 
input 

5.4.1. General considerations 

There are different ingredients required to calculate the HFM cross section with 
the formulas given in the preceding sections. Which ingredients impact different 
parts of the calculation in what manner is discussed below, but how a change 
in the transmission coefficients (of certain or all included transitions) affects the 
resulting cross section can only be understood with the help of ([T4l) . Similar to the 
determination of the energy-dependence of the cross section (which is crucial in the 
derivation of the relevant energy windows in Sec. [3]), the sensitivity of the cross 
section and rate to a change in the nuclear properties of the participating nuclei 
depends on which transmission coefficient (or averaged width) actually affects the 
cross section while the others cancel from ([M]) . As already pointed out in Sec.[3l the 
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discussion applies to both the BWF and the HFM. Thus, the sensitivities are rather 
well known in certain parts of the experimental community studying resonances and 
in the field of nuclear data evaluation. Since they do not seem to be so well known in 
Nuclear Astrophysics and their implications also have to be interpreted in terms of 
stellar cross sections and astrophysical reaction rates, it is helpful to outline briefly 
the main points here. 

By comparison to (p^ and (1551) we find that (|14l) transforms to 

rp^ rpFh rpAarpFb 

7^lab = and 7^,.atc = -^^-^ (67) 

^ Jit ^ Jit 

for each J, 7r-dependent summand. In laboratory cross sections (usually with /i = 0) 
only few J, 7r-summands contribute due to the spin selection rules, and so the appli- 
cation of 7?.iab to determine the sensitivities to changes in the widths is straightfor- 
ward. The situation is different for stellar rates and ratios 7?.ratc because transitions 
from excited target states additionally contribute and more terms in the sum may 
be relevant, depending on the nucleus and the plasma temperature. It is interesting 
to note that Tj°* includes in both cases. Depending on whether the average 
entrance transmission Tj^ significantly defines the size of T^, a variation of the 
entrance transmission will affect the total transmission more or less. This can lead 
to a different sensitivity of laboratory cross sections to the entrance channel than 
to the transmission in the exit channel, even if everything else is comparable. The 
interpretation of experimental results concerning the impact on stellar rates has 
also to proceed carefully. For example, if a strong dependence is found on the en- 
trance transmission coefficient and some deficiencies when comparing the model to 
experimental data, this does not necessarily mean that this is of relevance to the 
astrophysical rateO Since the astrophysical rate involves Tj^ also in the numerator, 
it may more often cancel with the denominator, even if it would not for Tj^ in the 
numerator. Only for rates with low SEF TZuh ~ 7?,rate- This shows again that stellar 
rates and laboratory cross sections do not have a one-to-one correspondence and 
additional, mostly theoretical considerations have to be included. 

Several special cases can appear in (1571) : (i) The larger of the two transmissions 
in the denominator also dominates the numerator; then the cross section or rate will 
change similarly to a change of the smaller transmission coefficient in the numerator 
and be oblivious to any others; (ii) Three channels will be important when neither 
of the transmission coefficients in the numerator dominates the total transmission; 
then any change in the two transmissions in the numerator will be translated into 
an equal change in the cross section but a change in the one determining the 
denominator will result in an inverse proportional change in the cross section. In 

^In this discussion, the term "transmission coefficient" can also be replaced by "width" which 
may be more familiar to some readers. However, if averaged widths (f) = T/{2iTp) or strength 
functions Sf = T/{2tt) can be determined experimentally, it is more useful to obtain the latter 
because they are directly proporti onal to the transmission coefficients T without an additional 
dependence on the compound NLD.E2] 
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Fig. 12. Sensitivities s for ®®Ru(p,7)^'^Rh when multiplying the transmission coefficients (averaged 
widths) by a factor of 2. 

the general case when T}°* is not dominated by a single channel the situation is 
more complicated as any change in a transmission coefficient will not fully affect 
the cross section at a similar level. 

A helpful visual aid to estimate the relative importance of the different channels 
is a sensitivity plot. I define the sensitivity s as a measure of a change in the cross 
section /o- — (Tncw/coid as the result of a change in a transmission coefficient (or 
width) by the factor fi^, with s = when no change occurs and s = 1 when the 
cross section changes by the same factor as the transmission coefficient (or width): 



Plotting s as a function of the cm. energy yields a plot like the example shown in 
Fig. [12] for the reaction ^^Ru(p,7)^''Rh. Its astrophysically relevant energy window 
is 1.63 < E < 3.42 MeV for the typical p-process temperature T = 2.5 GK.ISl It 
can clearly be seen in Fig. [T^lthat the sensitivities are very different at lower and 
higher energies. For example, a measurement closely below the neutron threshold 
would be in a region where s is largest for the 7 transmission coefficient (or width) 
but smallest for the proton transmission, just the opposite of what is found in the 
astrophysically relevant energy region. Above the neutron threshold the situation is 
even more complicated because there is additional sensitivity to the neutron channel 
(dominating r*°'), although not as large as to T ^^^^'>'. If any discrepancy between 
measured and predicted cross sections was found, it would be hard to disentangle 
the different contributions. In this example, no information on the astrophysically 
important proton transmission coefficient can be extracted from a measurement at 
higher energies. 

The above example also shows that some of the standard assumptions usually 



s = 




if > 1, > 1 or < 1, < 1 , 
if < 1, > 1 or /<x > 1, /a; < 1 • 



(68) 
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used in nuclear physics experiments do not apply. For instance, it is usually assumed 
that the 7— width is always smaller than particle widths and therefore a capture 
reaction will always only be sensitive to the 7— width. Many, if not most, astro- 
physical reactions with charged particles, however, proceed close to or below the 
Coulomb barrier. This leads to very narrow charged particle widths and they may 
well become smaller than the 7— widths. Even neutron widths close to a neutron 
threshold become smaller than photon widths (which are relatively independent of 
energy compared to the particle widths). Therefore it is always important to closely 
inspect the widths and to perform a thorough sensitivity study when investigating 
the astrophysical impacts of changes in the HFM inputs. 

From the general considerations above it also follows that it is advantageous 
to use reactions with neutrons in one channel for investigating the sensitivity to 
the charged-particle optical potential, i.e., using (a,n), (n,Q!), (p,n), or (n,p) reac- 
tions. Except within a few keV above the neutron threshold, the neutron width 
will always be much larger than the charged-particle width, even at higher than 
astrophysical energies. Therefore it will cancel with the denominator in Eq. (|67|) 
and leave the pure energy dependence of the charged-particle width. On the other 
hand, the neutron channel may not be open at the energies required to study the 
astrophysically important charged particle width and extrapolations have to be 
performed. Furthermore, this shows that it is difficult to obtain information on 
the neutron potential from reactions. But this also implies that the sensitivity of 
astrophysical rates to the neutron optical potential is not high. 

5.4.2. Relevance of nuclear input 

Nuclear properties and how they affect the calculation of transmission coefficients 
(or averaged widths) are briefly discussed in the following. This is by no means 
meant to be an exhaustive listing and discussion but rather some exemplary points 
are taken to explain the special requirements in the calculation of astrophysical 
rates and to point out the challenges. Further information on the nuclear input 
used in statistical model calculations for nuclear physics and nuclear astrophysics 
can be found, e.g., in Refs.|4l[9l[l3l[83l[26l[33l[9ll[9g[Sg[T03 references 
therein. 

While the HFM has been used extensively to study reaction data, its application 
in Nuclear Astrophysics has a slightly different focus. Standard nuclear physics 
investigations use the HFM by including measured or known properties of nuclei. 
Then reaction cross sections can be reproduced with high accuracy. Stellar cross 
sections not only require the modification of the HFM shown above but also the 
inclusion of mostly unknown further input because the reactions proceed at much 
lower energy and/or involve unstable nuclei. The required properties have to be 
predicted globally for a large number of nuclei and it has to be realized which 
properties are important in which astrophysical process. This is where the challenge 
lies and where basic research is necessary, going beyond a mere application of a 
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seasoned model. Although intellectually the most satisfying, modern microscopic 
models still are not applicable for large-scale predictions and do not produce all 
the necessary data with sufficient reliability. This calls for clever combinations of 
microscopic models and parameterizations (which can also use dependences derived 
from microscopic calculations for a limited range of nuclei) . Furthermore, additional 
information from experiments is required to constrain such models. 

5.4.2.1. Masses: Nuclear masses or, rather, mass differences determine separa- 
tion energies and reaction Q-values. In this way, they also determine the range of 
transition energies to be considered in each reaction channel and through this the 
relative importance of a channel. Except in the NSE equations p4p , masses always 
appear in mass differences. This poses a potential difficulty when a certain mass re- 
gion is not fully explored experimentally yet. Care has to be taken to avoid artifical 
breaks and structures in the mass surface when calculating mass differences from 
a mix of experimental and theoretical masses. This is usually considered in codes 
especially written for astrophysical applications. On the other hand, it is expected 
that mass differences can be measured as well as predicted with higher accuracy 
than single masses. This seems reassuring for the calculation of astrophysical rates 
far from stability. 

A change in the mass of a nucleus impacts the rates in two different manners. 
First, the separation energies are altered in the reaction channels including the 
nucleus. This leads to a change in the transition energies in these channels (see 
Fig. [2]). If the change is large, also more or fewer transitions may become possible. 
Although charged-particle transitions sensitively depend on the interaction energy, 
it has to be realized that even in this case the change in Q-value plus ejectile energy 
has to be considerable to have a sizeable impact on the rate. The entrance channel 
is not affected. 

Secondly, a change in the Q-value changes the relation between forward and 
reverse rate as shown in (f27| and (f28|) . It is highly sensitive to a small change in 
Q- value due to the exponential dependence. This impacts the temperature at which 
the two rates become comparable and at which equilibria are reached. On the other 
hand, if forward and reverse rates are different by many orders of magnitude (i.e. 
for large |Q|), there may not be much astrophysical impact, nevertheless. 

5.4.2.2. Properties of ground and excited states: It is standard procedure 
to include spin, parity, and excitation energy of low-lying discrete levels when cal- 
culating the transmission coefficients. Information about discrete states comes from 
experiments or from nuclear theory (single particle states and shell model states). 
Close to stability, a large number of excited states are well known. However, often 
excited states of nuclides not produced in reactions on stable target nuclei are only 
partially known experimentally, even if the nuclide itself is stable. Obviously, the 
situation worsens further off stability. Therefore, nuclear spectroscopy is important 
to provide the database for reaction modelling. Discrete states are not only impor- 
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tant in the HFM but also for resonant reactions treated in an i?-matrix or BWF 
approach (unbound states) and even more for the treatment of direct reactions (see 
Sec. 15. 5|) (bound states). 

It is very important to have a complete set of excited states because a large 
number of missing levels at and below a given excitation energy would lead to an 
incorrect prediction. Therefore it is important to set a cut-off in the excitation 
energy above which no states are included (even if there are some data) and below 
which the assumption of a complete level set holds. Above the cut-off, a NLD 
is employed (see below). It is difhcult to define a useful cut-off energy but blindly 
including all existing data leads to worse results than neglecting too many, provided 
a reliable NLD description is used. Usually, an educated guess has to be made by 
comparing the level data to NLD predictions. 

5.4.2.3. Nuclear level density: It should be noted that the level density p used 
throughout this paper follows the "experimental" definition of number of observed 
levels per energy interval p{E) — j tt t J i ^ where p{E,J,Tr) is the observed 
number of levels with spin J and parity tt in a small energy interval around an 
excitation energy E. It is not to be confused with a state density uj{E) = X)j7r(2>^ + 
l)p{E, J, tt) appearing in microscopic nuclear theory, such as the shell model. The 
terms "level density" and "state density" are used inconsistently in the literature. 
For the relation of the two types of density, see, e.g., Refs. 11031 1771 

In the regular HFM, the nuclear level density only enters in the calculation 
of the transmission coefficients (Eq. I64p when there are no or not enough discrete 
states known and therefore its importance depends on how many low-lying, discrete 
states were included (see above). 

Because of the energetics connected to particle emission (see Fig. [2]), usually 
only a small fraction of the particle width is due to transitions calculated with 
a level density at astrophysically relevant projectile energies as long as discrete 
excited states are known. This is especially true for neutron-capture reactions due 
to their low interaction energies. Reactions with charged particles in entrance or 
exit channel prefer somewhat higher projectile energies due to the location of the 
relevant energy window (Sec. |3] and thus may show slightly larger sensitivity to the 
level density in the target and final nucleus. (These sensitivities are different from 
those encountered in usual nuclear reaction studies where projectile energies are 
much larger.) 

On the other hand, there is a larger range of energies available for 7-transitions 
(see Fig. [2] grey areas signify transitions calculated by integration over a level 
density, as shown in the second term on the right hand side of Eq. I64p and thus the 
impact of the level density in the compound nucleus will be largest. In most cases, 
it is accurate to assume that a variation of the level density will mostly (or only) 
affect the 7- widths. 

Generally, the impact of a change in the NLD will be larger for nuclides with 
fewer known discrete states, i.e. far from stability. On the other hand, the Q- values 
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of the astrophysically relevant reactions become lower and thus ground state tran- 
sitions to the final nucleus may again dominate. 

It is instructive to remember which excitation energies are the most important 
ones. For particle transmission coeflBcients the transitions with highest relative en- 
ergy are most important, i.e. to the ground state and the lowest excitation energies 
in both target and final nucleus (similar considerations apply to direct reactions). 
Discrete levels or a NLD have to be known there. For electromagnetic transitions, 
the relevant excitation energy is around i?form/2 in the compound nucleus (see the 
discussion of electromagnetic transmission coefficients below and Ref. 11041 for a 
detailed explanation). 

Due to the low excitation energies involved it is important to a) include the cor- 
rect ground state spin and parity, and b) to account for a possible parity dependence 
of the NLD at low excitation. Both is automatically ensured when experimental in- 
formation is available up to sufficiently high excitation energies. Otherwise, ground 
state properties and NLDs have to be predicted. This introduces additional uncer- 
tainties in rates far from stability. 

Until recently, astrophysical rate predictions made use of equally distributed 
parities p+ = p- = Modern rate predictions include parity-dependent 

NLDs in yarious man ners. Either, microscopic NLDs are directly ut ilized in 
calculationPSHniMI or a parity-dependence is applied to a total NLD.I^^ ^^ I ^^ I 
The advantage of the latter approach is that it can conveniently be applied to total 
NLDs from different sources and for a large number of nuclei. The total NLD p 
is not changed but the parities are redistributed according to excitation energy. 
This approach was used to study the impact of a parity-dependence across the nu- 
clear chart. Since capture reactions mainly populate higher lying states, for which 
an equipartition of parities already is a good assumption, the impact of a parity- 
dependent level density is small, unless very low inform is encountered due to low 
projectile separation energies and low plasma temperatures.!^ According to the 
discussion above, however, the impact is larger in particle emission channels and 
for direct reactions. 

A modified HFM is introduced in Sec. 15.3.31 above, which accounts for the rela- 
tive level distribution at the compound formation energy -Eform- Again, the impact 
of a parity-distribution is small for sufficiently large inform and only beco mes impor- 
tant far from stability at low separation energies (in the r-process).!^ However, it 
remains doubtful whether this actually is of astrophysical relevance as it is not clear 
whether the involved nuclei can be produced outside of equilibria (see Sec. l4.2|) and 
whether the uncertainty introduced by using the statistical model for nuclei with 
such low NLD at inform is not much larger than the impact of the parity-dependence. 
On the other hand, a dependence on the total NLD and on the J distribution in 
the further modified HFM may be more important at a larger range of excitation 
energies and thus also for nuclei closer to stability (see Sec. I5.3.3| . 

In any case, the uncertainties introduced by the predicted NLDs in rates far 
from stability are overall much smaller than the uncertainties stemming from other 
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input to the HFM calculations, such as optical potentials and photon transmission 
coefficients. 

5.4.2.4. Optical potentials: Optical potentials are required in the solution of 
the radial Schrodinger equations to determine the particle transmission coefficients 
as shown in (j65l) . Together with the electromagnetic transitions strengths (see be- 
low), the unknown optical potentials at low energy give rise to the largest uncer- 
tainties in astrophysical reaction rate predictions. There is a combination of two 
problems involved in the determination of appropriate optical potentials for astro- 
physical applications: the prediction of optical potentials for highly unstable nuclei 
and the extension to the astrophysically relevant energies. There is a large amount 
of reaction and scattering data along the line of stability and many parameteriza- 
tions (usually of the Saxon- Woods shape), partly mass- and/or energy-dependent, 
are available. However, most scattering experiments to derive optical potential pa- 
rameters have been performed at several tens of MeV, far above the astrophysically 
relevant energy window (see Sec. [3]). Even at stability, there are almost no data (not 
even reaction data) for charged-particle reactions at astrophysical energies. Measur- 
ing low-energy cross sections for charged-particle reactions is especially problematic 
due to the Coulomb barrier causing the astrophysically relevant cross sections to be 
tiny. Even more problematic is the standard way to obtain information on optical 
potentials through elastic scattering experiments. The scattering cross section at 
low energy becomes indistinguishable from Rutherford scattering. 

In the optical model, an imaginary part of the potential appears whenever there 
is loss of flux from the elastic channel due to any kind of inelastic process. The 
reaction cross section is sensitive to both real and imaginary part of the optical 
potential because they determine the relation of real an d imaginary pa rt of the nu- 
clear wavef unction and thus the transmission coefficient! ^^ ! '*^^ ! '*^^^ ! ^ ^ Microscopic 
approaches to derive the optical potential are preferrable over parameterizations 
when predicting rates far from stability, especially because the available parame- 
terizations were derived at far too high projectile energies. Especially the imaginary 
part of the potential may vary strongly with energy, due to the energy-dependence 
in the available reaction channels included in the imaginary part. Nevertheless, any 
more sophisticated theoretical approach also includes some parameters which have 
to be constrained by comparison with experiment and so even in this case there 
may be uncertainties at low energy (see below and the further discussion of how 
additional reaction mechanisms impact the optical potential in Sec. 15.5.1]) . 

An optical potential widely used for interactions of nuclei with neutrons and 
protons uses the Briickner-Hartree-Fock approximation with Reid's hard core 
nucleon-nucleon interaction and adopts a local density approximation.'^JJ^ A low- 
energy modification of this potential was provided specially for astrophysical 
applications.^^ The latter has become the standard potential in predictions of 
astrophysical rates and is generally very successful when compared to the scarce 
experimental data at low energy. As stated in Sec. 15.4.11 astrophysical rates are 
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Fig. 13. Astrophysical S-factors for (p,7) reactions on ^"Ge (data are from Ref. IllST i com pared 
with theoretical values obtained with dif feren t optical potentials: "standard" potential (JLM)(li2J 
modified standard potential (mod JLM)|3fi! and the potential from Ref. 1 11 7l (Bauge) . 
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Fig. 14. Same as Fig.[T3]for ^*Se(p,7) 



rather insensitive to the neutron potential. Regarding rates involving charged par- 
ticles in entrance or exit channel, an important difference to cross sections at higher 
energies is in the fact that astrophysical cross sections are mostly sensitive to the 
charged particle widths instead of the 7- or neutron widths. A series of (p,7) and 
(p,n) reactions was measured close to astrophysical energies recently (see, e.g., 
Refs.l36 l 1113(11141 and references therein). The latter reactions are especially useful 
for testing the proton potential because the neutron width will (almost always) be 
larger than the proton width at all energies. 

Despite of the overall good agreement when using the standard potential, cer- 
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Fig. 15. Astrophysical S-factors of *^Rb(p,n)^^Sr (exp. data from l36l l compared with theory using 
different proton pot entials and nuclear densities. Shown are results with nucle ar d ensity from 
a droplet modellSSl and from a HFB mod el w ith Skyrme interaction (HFB-02j -^^-^l, applied in 
the calculation of the "standard" potentialDS] (JLM) and a modified version of this potential 
(modJLM) with increased absorption. 
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Fig. 16. Astrophysical S-factors of *^Rb(p,n)*^Sr (exp. data froml^Sl) compared with theory using 
different prot on p otentials and nuclear densities. Shown are results with nuclear density from a 
droplet modeffSl] ^nd from a HFB model with Skyrme interaction (HFB-02 j^^^ l, applied in the 
calculation of the potential of Ref. 11161 (Bauge) and Ref. '117' (BaugeLane). 



tain systematic deviations at low energy were found recently (see, e.g., Figs. fTSlfMl 
more examples are shown in Ref. lll3p . It was found that an increase in the strength 
of the imaginary part at low energies considerably i mproves th e reproduction of the 
data (denoted by "mod JLM" in Figs. [T3 l [M l [T5| ! ^" l ^^ tll5| increased absorp- 
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tion is permitted within the previous parameterization because the isoscalar and 
especially the i sovector compon ent of the imaginary part is not well constrained 
at low energies! ^ ^^ I -*^-*^^!-*^-*^^ ! ^^^ ! Therefore, the change has to be energy-dependent 
(i.e. acting only at low energy). More low-energy data is required to obtain a better 
picture. 

Figures [T3l [Ml also shows results obtain ed w ith another recent. Lane-consistent 
new parameterization of the JLM potential.'!!-^' Although it showed improved per- 
formance at higher energy, it yields worse agreement at astrophysically low energy. 
This is understandable as neither does it include the additional modifications of 
Ref. I112( nor ca n it coiis train well the low energy part because it was fitted to data 
at higher energy.lH^HI^ Simil ar considerations apply to another recent reevaluation 
of the standard potential.l^l^ 

Required input to the calculation of this type of optical potentials is the nu- 
clear density distribution p™'^**. Figure ITSl shows results when employing a dro plet 
model density -'^^'^ and one from a Hartree-Fock-Bogolyubov model (HFB-2).1SU 
For the reactions considered here, the droplet description yields better agreement 
to the data in both absolute scale and energy dependence of the theoretical S- 
factor. For comparison. Fig. [16] also shows the results when employing the optical 
potentials of R efs. II 161 1 1 1 71 with both densities. In the original work, HFB densities 



were employed! 1 1 ^ 1 H ZI pgj- further information on how nuclear density distributions 
affect the astrophysical rates see the subsection on density distributions below, in- 
cluding Figs. [111-1271 

The standard potential for neutrons and protons seems to work very well com- 
pared to the situation encountered when exploring the adequacy of a-|-nucleus op- 
tical potentials for astrophysics. Global parameterizations describing scattering, re- 
action, and decay data have been notoriously hard to find for a-particle potentials. 
Somewhat surprisingly, a mass- and energy-independent potential of Saxon- Woods 
type has been quite successful an d is widely used to evaluate reaction data and 
also for astrophysical applications The potential was fitted to scattering data 
at 26.7 MeV a-energy for a wide range of nuclei. However, it became obvious early 
on that it may be impossible to find a global potential with a pre dictive po wer 
comparable to those for nucleonic projectiles, especially at low energy.E^HlQl xhe 
number of optical p otent ial parameters can be reduced by using folding potentials 
Uf for the real part,^!^ 



C/F(r) = A j d'rA j d'r,pr'\rA)pT\r,)x 

v,B {E, p'"^" = pr" + pr'\S = \r + fa - fAl) . (69) 

In this expression r is the separation of the centers of mass of the two interacting 
nuclei, p™^" and p™*^** are their respective nucleon densities and A is an adjustable 
strength factor. The factor A may differ slightly from Unity because it accounts 
for the effects of antisymmetrization and the Pauli principle. The effective nucl eon- 
nucleon interaction v^s for the folding procedure is usually of the DDM3Y type.^^^ 
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Fig. 17. Astrophysical S-factors of l'''*Sm(a,7)i''*Gd (exp. data from Ref.[127j) compared to HFM 
calculations with different optical a potentials, from Refs. [1271 (fit). fllQI (McFadden*). 11381 (Froh). 
11321 (Avri 2010). Transmission coefficients were computed with the routine by Refs. [1461 11471 



Density distributions have to be taken from experiment or theory (see the section on 
matter density distributions below) . A global parameterization of the real p art with 
such folding potentials was found based on extensive scattering data.L!^^ Unfortu- 
nately, there is no simple description for the imaginary part, for which shape and 
strength have to be energy-dependent. Correlations with the compound NLD, dif- 
ferent parametrizations for the energy-dependence of the strength (Fermi-function, 
Brown- Rho dependence) , and an energy-dependence in the relative strengt h of vol- 
ume and surface imaginary parts have been suggested.^-^^ -'^^^ ^"^^ -'^27 130 131|1 32] 
Again, those extrapolations to low energy are only loosely constrained due to the 
lack of scattering data. Based on apparently different potentials required for the 
description of a-particles in entrance and exit channel and the fact that a poten- 
tial fitted to reaction data is able to describe a number of reactions but does not 
reproduce scattering data, it was suggested that there may be some dependence on 
nuclear temperature and that absorption and emission potentials may be different 
from scattering potentials. "'^'^"'^ ^^'"^ ^'^'^ ^"^^ Moreover, at energies close to or below 
the Coulomb barrier, so-called "threshold anomalies" have been observed, a rapid 
variation of optical potential parameters with energy.tlSSl It has been shown that 
the dispersion relation conne cting real and imaginary part of the optical potential 
is essential to describe these.l^^^ There is a large literature on different local and 
global parameterizations of a-|-nucleus potentials, mostly at high energy, underlin- 
ing the lack of a coherent nuclear physics treatment. A complete review cannot 
be provided here but see, e.g., Refs. [71 |98l ESI I120[ 11311 and references therein for 
further details. 

For astrophysical application s, early rea ction rate predictions made use of an 
Equivalent Square Well potential! I ^"^ ! -*^"^^ Later on, the potential of Ref. 11191 was 
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Fig. 18. Same as Fig. 1171 but transmission coefficients were computed with a modern routine for 
the solution of the Schrodinger equation at energies close to the Coulomb barrier (Refs. [260ll93ll . 

used. With accumulating evidence that low-energy data deviates from the predicted 
cross sections more complicated parameterizations were tried, based on various 
combinations of scattering and reaction data, but these did not lead to a consis- 
tent picture. On the other hand, the potential by Ref. IT?81 (see also Refs. IMl 11391 
I140p was fitted to simultaneously reproduce low-energy reaction cross sections of 



Nd(n,a)i''°Ce, ^^'^Sm{n,ay^'^m, and "4Sm(a,7)"8Gd. Although the potential 



does not describe scattering data (at higher energy) it was found to work surpris- 
ingly well for low-energy cross sections for target nuclei across a large mass range 



potential of Ref. 11191 it predicts systematically lower cross sections at low energy 
due to a shallower real part. Although not fully satisfactory yet, this shows that the 
main change required by a global a-potential is to reduce the predicted low-energy 
cross sections by factors of about 2 — 3. 

The only and remarkable exception to these factors known so far is in the com- 
parison to experimental data for ^^^Sm(a,7)^^®Gd.lJ^^This reaction is im portant in 
the astrophysical context not just generally in p-process calculationJl^ but specif- 
ically in deriving astronomical timescales from the Nd/Sm abundance ratios mea- 
sured in meteoritic inclusions or, vice versa, to determine production ratios of these 



10.2 — 13.0 MeV, the astrophysically relevant energy window extends from 9 MeV 
downwards. Above 11.5 MeV the S-factors obtained with the potential of Ref. [TT9l 
are too high by a factor of three but the data is well described wi th th e potential of 
Ref. 11381 and a potential derived from scattering data at 20 MeV.^^^^ However, be- 
low 11.5 MeV the energy-dependence changes and requires further modifications of 
the potentials. At the lowest measured energy of 10.2 MeV, the measured S-factor 
differs by a factor of 10 from the one predicted with the potential of Ref. '1191 In 




elements 




astrophysical S-factors were measured from 
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Fig. 19. Dependence of predicted astrophysical S-factors of ^''^Sm(a,7)^''*Gd on the Coulomb 
radius parameter used with the potential of Ref. 11271 

Ref. 11271 a potential with energy-dependent imaginary part was fitted to describe 
the data. It predicts S-factors at and below 9 MeV which are orders of magnitude 
lower than those predicted with the global potentials. Although the original problem 
was a too low calculated Nd/Sm ratio as compared to what is found in meteorites, 
such low S-factors yield much too high ratios and add itional astrophysical dilution 
effects have to be invoked in an ad hoc manner .ll^^ The S-factors are shown in 
Fig. [171 where also new results obtained with the recent potential of Ref. 1321 are 
included. It should be noted that this reaction may be a special case because it 
is strongly endothermic. Nevertheless, it was shown th at t he stellar enhancement 
factor is lower in the capture direction (see also Sec. l4.3p .l^ Since the extrapolation 
to astrophysical energies strongly hinges on the data points at the lowest energies, 
an independent remeasurement of this reaction at comparable or lower energies is 
highly desireable. 

Further issues in the determination of S-factors below the Coulomb barrier are 
illustrated in Figs. [TTl - [20l Figure [17] shows the results obtained when using the 
routine of Ref. [T46| (as also used in Refs. [Ml U UMl [HIl [HZl [US]) for the solution 
of the Schrodinger equation and the determination of the a-particle transmission 
coefficients. It does not work properly when applied at energies at or below the 
Coulomb barrier as a comparison to results obtained with a modern routine shows 
(Fig, lisp. The discrepancy between predictions and measurement is even enhanced 
when using an appropriate method. Since the S-factors are very sensitive also to the 
wavefunction far outside the nuclear center, they show a strong dependence to the 
shape and width of the effective barrier (determined by the sum of nuclear potential 
and Coulomb potential). Results obtained with the optical potentials of Refs. 11271 
and [1321 with fine-tuned energy-dependences of the real and/or imaginary part also 
have a strong sensitivity to the Coulomb radius parameter, as shown in Figs. llQl and 
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Also, when using the potentials of Refs. 11191 11381 there is no sensitivity to the 
Coulomb radius in the investigated energy range. These additional comphcations 
show that any extrapolation to subCoulomb energies has to be performed very 
carefully and that it is difficult to construct a global potential. 

Further experimental data (on scattering and reactions) are especially in de- 
mand for improving the optical potentials (see also the discussion of additional re- 
action mechanisms contributing to the absorptive part of the potential in Sec. l5.5.T|) . 
Currently, progress is hampered by the lack of systematic reaction (and scattering) 
data at astrophysically relevant energies, even at stability. 



5.4.2.5. Electromagnetic transitions: First, a few words on the energies of 
the emitted 7-rays and their significance for changes in 7-ray strength functions 
(drawn from experiment or from theory) are in order. As shown in Fig. [21 the en- 
ergies of emitted 7-rays are in the range < < E'form- Therefore, the behavior 
of 7— strength functions at low energy have to be known. Since the strength of 
the 7-transition scales with some power of Ej^ 7-transitions with higher energies 
are favored. On the other hand, the number of available endpoints of the transi- 
tions increases with increasing excitation energy of the nucleus because the NLD 
increases rapidly. This competition between transition strength and NLD gives rise 
to a peak in the 7-emission energies as shown in Fig. 1211 This peak is fragmented 
when certain transitions to discrete excited states are dominating. This is mainly 
the case far off stabi lity f or captures with low Q-values, forming a compound nu- 
cleus with low NLD.1^ Figures [22l [23] show examples for the 7-energies which 
maximally contribute to the reaction rate integral. Interestingly, it can be seen that 
for astrophysically relevant projectile energies, the 7-energies with the strongest 
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impact are between 2 — 4 MeV unless the level density i s so low that only few 
transitions are allowed (usually for highly unstable nuclei) Then the relevant 

is = i^proj + ^-scp — E'^ (in the figures this is the ground state with E^ ~ 0) 
instead of the almost constant value below inform- 

The 7-emission peak defines the range of 7-energies at which changes in 
the strength function have largest impact as well as the excitation energies at 
which the NLD is most important. This also holds for the reverse reaction (pho- 
todisintegration) under stellar conditions because the additional, linear weight 
1 — E'^/{Eioim ~ E^) from p2)) has a much weaker E^ dependence then both the 
NLD and the 7-strength. In consequence, changes in the strength function around 
this energy have the largest impact. Testing strength function models outside the 
energy range will not be relevant to astrophysics. Unfortunately, such low energies 
cannot be probed by photodisintegration experiments because they are below the 
particle separation energies, at least close to stability. Such experiments would al- 
low to study strength functions in the most direct way (but they cannot test the 
astrophysically relevant rates, either, see Sec. 14. 3p . Other types of experiments are 
complicated by the fact that the observables are generated by a convolution of 
different nuclear properties (such as the dependence on the NLD or different spin 
selectivities of transitions) which have to be known and disentangled. 

At least the dominant 7-transitions (El and Ml) have to be included in the cal- 
culation of the total photon width for astrophysics. Some codes offer the possibility 
of including higher order transitions. There are two issues involved: 1) Obtain- 
ing, understanding, and modeling photon strength functions (PSF) at stability; 
2) predicting strength functions far from stability by using parameterized or mi- 
croscopic models, predicting the nuclear properties (e.g., deformation) entering the 
descriptions. Despite of decades of experience in studying nuclear reaction data, the 
understanding of the electromagnetic transitions between nuclear states i s lirnited 
and predictions are subject to considerable uncertainties, even at stabilityl ^^^ l ^^^ l 
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The paradigm in the field of studying electromagnetic transitions is the validity 
of the reciprocity theorem (jl7p also when applied to photon emission and absorp- 
tion, and the independence of the PSFs from the nuclear structure of the initial 
and final stat es (e xcept for spin and parity, selecting the allowed multipolarity of 
the radiation) ffffl This is called the Brink hypothesis! -*-^-*-! This is also the basis for 
the construction of the effective cross section (IT9|) . the expression for the reverse 
rate and the introduction of equilibrium abundances (Sec. l4.2| ). The Brink hy- 
pothesis has been studied extensively in experiments and its violation would have 
grave consequences not only in nuclear reaction theory but also for astrophysical 
reaction rates and network calculations. 

Among the collective modes of nuclei the electric dipole (El) excitation has the 
special property that most of its strength is concentrated in the isovector giant 
dipole resonance (GDR). Macroscopically, this strong resonance is described as a 
vibration of the charged (proton) matter in the nucleus against the neutral matter 
(neutrons) . The transmission coefficient in a nucleus with charge number Z, neutron 
number iV, and mass number A — Z + N can be parameterized aJ^^ 

Tbi(£;^) = -— 2 V ---2 — ^2 ""xT , ^T^2 ^ ■ (70) 



8NZ l + X ^ i rGDR,,;-E; 

3 ^ rn^^2 A. 3 (i52 _ ^2 ^)2 + r^oR,.^^^ 



Here, uip is the p roton mass, x(= 0.2) accounts for the neutron-proton exchange 
contribution,H^^and the summation over i includes two terms which correspond to 
the split of the GDR in statically deformed nuclei, with oscillations along {i = 1) 
and perpendicular {i — 2) to the axis of rotational symmetry. In this deformed 
case, th e two resonance energies are related to the mean value calculated by the 
relationJiSlI 

EgDRS + 2i?GDR,2 = 3i?GDR , 

i^GDR,2/-BGDR,i = 0.91177 + 0.089 . (71) 

The deformation parameter r/ is the ratio of the diameter along the nuclear symme- 
try axis to the diameter perpendicular to it, and is obtained from the experimentally 
known deformation or from mass model predictions. Many microscopic and macro- 
scopic models have been devoted to the calculation of the GDR energies -Egdr and 
widths Fgdr- 

Of special interest here is the low-energy tail of the GDR. It is a long-standing 
question of nuclear physics to specify how much of the El strength is still present at 
energies far below the GDR maximum, which also encompasses the astrophysically 
relevant energy region. Theoretically it has been shown that it is justifie d to describe 
the GDR by a Lorentzian also below the particle emission thresholds ) 100 | 150 | 155 | 
Various experimental attempts to determine the low-energy extension of the GDR 
for heavier nuclei have led to conflicting results. Neutron-capture experiments of- 
ten have indicated an overshoot of the Lorentzian over the observed El strength 
at the low-energy 

tailoftheGDR.™l56ll57ll58] 

On the basis of these data theo- 
retical explanations have been proposed to explain the differences ) 100 | 159 | 160 | 15 H1 
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Fig. 22. Maximally contributing 7-energies E-^ compared to the projectile separation energies i/scp 
when capturing i?proj = 60 keV neutrons on Sn isotopes. The mass number A is the one of the 
final (compound) nucleus. See Ref. 11041 for details. 




Fig. 23. Same as Fig. [22] but for Bp^oj = 5 MeV protons. See Ref. lTOil for details. 



Photon-scattering exp er imen t s , h-owever , are in m any cases in good agreement with 
the Lorentzian extrapolation.^^ ^ l ^^^ l ^^^ l ^^^ l ^^^ l Unfortunately, they do not access 
the astrophysically relevant energ ies. A lso other experiments are used to extract 
PSFs, e.g., •^He-induced reactions.^^ Recent investigations have sho wn that the 
El strength can be described by Lorentzians for a large range of nuclei .^^^^^^^^^^^ 
It was also shown that indirect determinations of PDFs are prone to larg e uncer- 
tainties due to the experimental difficulties and s ome claims of enhancementP^^^l^ 
at low have been premature.'lSQllSHIIMIMl 

However, in some nuclei extra strength at low energy with respect to the 
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smooth Lor entzian was found consistently and denoted as pygmy dipole resonance 
(PDR) IlSSl The PDR has experimentahy been studied so far in sphe rical nuclides 
around Z,N = 20, 28, Z = 50,N = 82, and in the doubly magic ^ospb.H^ Theoret- 
ical approaches describe the PDR as caused by an oscillation of excessive neutrons 
against the symmetric proton-neutron system (see, e.g., Refs. 11711 1172[ 11731 1174|) . 
Other oscillation modes were also proposed (e.g., scissor modes) which may also 
add strength beyond the Lorentzian tail of the GDR. In any case, extra strength 
in the low-energy tail only has an astrophysical impact if it is within the relevant 
energy range defined above. Although a PDR may le ad to an increase by several 
orders of magnitude in the astrophysical capture rate,^I^3IIl5] ^j^jg depends sensi- 
tively on its location and width. Different models give varying predictions of these 
crucial properties. Depending on the microscopic model used, the pygmy resonance 
is sometimes pr edicted at too hig h or too low an energy as to have any astrophysi- 
cal consequence! 1 1 77 | 1 78 | 1 79j investigation of this issue is required. It is 
important to note that the uncertainties in predicting the PDR enter additionally 
to the general uncertainties still present in the prediction of GDR energies and 
widths. Together with the predictions of optical potentials for the particle trans- 
mission coefficients, these are among the largest uncertainties in the determination 
of astrophysical rates. 

PSFs of higher multipole order are even less studied due to their small contribu- 
tions to cross sections. There are several descriptions available for Ml transitions, 
starting from PSFs from the simple single particle approach to more sophisticated 
(but less thoroughly tested) models! I "^ l ^^ l "^" ! Unless energy-independent PSFs 
are employed, a relevant energy window similar to the one for El transitions will 
arise. 



5.4.2.6. Isospin: Isospin conservation restricts transitions to certain final states 
with the same isospin as the initial and compound states, i.e. AI — 0. Isospin 
conservation is not absolute and cross section measurements of isospin-forbidden 
reactions give an estimate of the size of the isospin breaking (or isospin mixing). 
Internal isospin mixing due to the Coulomb interaction and external mixing via the 
other reaction channels have to be distinguished. The HEM equation as shown in 
(|62|) with the transmission coefficients (|64l) does not account for isospin conserva- 
tion unless it is included in the transmission coefficients. In other words, complete 
isospin mixing is assumed. The calculation of the transmission coefficients can be 
generalized to explicitly treat the contributions of the dense background states with 
isospin = I^-^- and the isobaric analog states with I> — I< + l.^'^ 186 187| 
In reality, compound nucleus states do not have unique isospin and for that reason 
an isospin mixing parameter J, was introduced which is the fraction of the width 
of states leading to transitions."'^^'' For complete isospin mixing fx^ — I, 
for pure states fj,^ = 0. In the case of overlapping resonances for each involved 
isospin, fi^ is directly related to the level densities — p and p^, respectively. 
Isolated resonances can also be included via their internal spreading width F"'' and 
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a bridging formula was derived to cover both regimes!!^ 

In order to determine the mixing parameter ^.^ = ^^{E), experi menta l infor- 
mation for excitation energies of levels can be used where av ailable .^^^ Exper- 
imental values for spreading widths are also tabulatedff^^^^^ Inspection of the 
tables shows that internal mixing dominates and that the associated spreading 
width is nearly independent of r nass nu mber and excitation energy, facilitating the 
extrapolation to unstable nuclei.E™ Similarly to the standard treatment for the 

states (the regular transmission coefficients as shown above), a NLD can be in- 
voked above the last experimentally known level. Since the states in a nucleus 
{Z,N) are part of multiplet, they can be approximated by the levels (and NLD) 
of the nucleus {Z—1,N+1), only shifted by a certain energy E^- This displacement 
energy can be calcul ated and it is dominated by the Coulomb displacement en- 
ergy Ed — -I- e{inni Thus, the uncertainties involved are the same as in the 
prediction of the NLD and discrete excited states. 

The inclusion of the explicit treatment of isospin has two major effects on sta.- 
tistical cross section calculations in astrophysics which will be discussed 

belowfSSl 

the suppression of 7-widths for reactions involving self-conjugate nuclei and the 
suppression of the neutron emission in proton-induced reactions. Non-statistical 
effects, i.e. the appearance of isobaric analog resonances, can be included in th e 
treatment of the mixing parameter /i^ but will not be further discussed here.^lSSl 

The isospin select ion rule for El transitions is A/ = 0,1 with transitions — > 
being forbidden !^^^^ An approximate suppression rule for AJ = transitions in 
self-conjugate nuclei can also be derived for Ml transitions.!!^ 

In the case of (01,7) reactions on targets with N = Z, the cross sections will 
be heavily suppressed because I — 1 states cannot be populated due to isospin 
conservation. A suppression will also be found for capture reactions leading into 
self-conjugate nuclei, although somewhat less pronounced because 1=1 states 
can be populated according to the isospin coupling coefficients. This cross sec- 
tion suppression can be implemented as a suppression of the photon transmission 
coefficient. Some older reaction rate calculations treated this suppression of the 
7-widths completely phenomenologically by dividing by rather arbitrary factors of 
5 and 2, for (01,7) reactions and nucleon capture reactions, respectively.l^^ESIinS] 
This can be improved by explicitly accounting for population and decay of 
and states, and considering isospin mixing by the parameter ^^.^^"^ An as- 
trophysically important reaction with pronounced isospin suppression effect is the 
reaction ^°Ca (a,7)'^^ Ti which is responsible for the production of the long-lived ''^Ti 
in supernovaeP^^QSi Decay 7-emission of ^^Ti is observed in supernova remnants 
and can be used to test supernova models.l^^ 

Furthermore, assuming incomplete isospin mixing, the strength of the neutron 
channel \ yill be supp ressed in comparison to the proton channel in proton-induced 
reactions! I ^ ' This leads to a smaller cross section for (p,n) reactions and 
an increase in the cross section of (p,7) reactions above the neutron threshold. 
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Fig. 24. Ratios of (p,7) rates at T = 3 GK on stable and proton-rich nuclei, obtained by 
using the HFB-02 and HFB-14 nuclear density distributions in the calculation of the optical 



as compared to calculations neglecting isospin (i.e. implicitly assuming complete 
isospin mixing with fi'^ = 1). The isospin mixing p aram eter was varied in the theo- 
retical investigation of a ^^V(p,7)^^Cr experiment .^^^ It was found that complete 
isospin mixing closely reproduced the measured cross sections when width fluctu- 
ation corrections were considered. Width fluctuation corrections affect the (p,7) 
cross sections above as well as below the neutron threshold, whereas incomplete 
isospin mixing only reduces the cross sections above the threshold. Thus, the two 
corrections can be discriminated. Mainly from this result it was concluded that - 
contrary to width fluctuation corrections - isospin can be neglected. However, a 
closer investigation of the levels in ^^Cr (using results from Refs. 11881 14| shows 
that isospin mixing should be rather complete already at the_ neutron threshold 
(since the first state is almost 1 MeV below the threshold). ^"^^ This is also true 
for lighter targets. For reactions on heavier nuclei {Z > 30), however, the neutron 
and proton threshold, respectively, will still be in a region of incomplete isospin 
mixing and therefore isospin effects should be detectable there. On the other hand, 
this effect is not as important in the calculation of astrophysical reaction rates as 
the suppression of the 7-width because of the averaging over an energy range in the 
calculation of the rate, washing out the cusp effect. 

5.4.2.7. Nuclear matter density distribution: The density distribution p™^" 
of neutrons and protons inside a nucleus is needed to calculate the optical potentials 
for some choices of potentials (see above). T here are charge density distributions 
available from electron scattering experiment J ^^^ ^ on stable nuclei but the majority 
of density distributions for application across the nuclear chart comes from theoret- 
ical predictions. Rates are mostly sensitive to optical potentials, and thus densities, 



potentials QSniZI 
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Fig. 25. Ratios of (p,7) rates at T = 3 GK on stable and proton-rich nuclei, obtained by using 
the HFB-14 and droplet model nuclear density distributions in the calculation of the optical 
potentials.tlSMSSj 



at large radii because of the low astrophysical energies. 

Density distributions are required in the determination of optical potentials and 
the potentials used in direct capture calculations. Some optical potentials rely on 
nuclear density distributions from a certain model which w ere used when fitting the 
remaining open parameters to experimental data.EHUll 

Although modern microscopic models have considerably improved in predicting 
nuclear masses and radii, the differences between different models still are large. 
Because of the sensitivity of the transmission coefficients to the optical potentials, 
one would expect that small differences in the nuclear density distributions can give 
rise to large differences in the rates, especially at low plasma temperatures. As an 
example, Figs. [24l[25l show comparisons of (p,7) rates obtained when using densities 
from Hartree-Fock-Bogolyubov (HFB) models and fro m the drople t model in the 
JLM potentials while leaving all other input unchanged ) 181 | 182 | 183 | Densities from 
HFB-2 were included in the Recommended Input Parameter Library for Hauser- 
Feshbach calculations RIPL-2 (Ref.'98|) and the IIFB-14 densities are in its successor 
RIPL-3 (Ref. [99|) . The left panel shows the comparison of rates obtained within the 
"family" of HFB models. The right panel shows a comparison of rates obtained with 
droplet model and HFB-14 densities. While the ratios stay well within a factor of 
two even with the seasoned droplet densities, it is interesting to note that they 
reach Unity when approaching the proton dripline. 

Figures [26l [27l are the same as above but for (n,7) rates at T = 1 GK with stable 
and neutron-rich target nuclei. Here, the largest ratios appear for very neutron-rich 
nuclei but again recede to Unity when approaching the neutron dripline. Overall, 
the maximal ratios are higher than for proton capture but stay within a factor of 
ten. This is partly due to the fact that nuclei further from stability are involved 
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Fig. 26. Ratios of (11,7) rates at T = 1 GK on stable and neutron-rieh nuclei, obtained by 
using the HFB-02 and HFB-14 nuclear density distributions in the calculation of the optical 
potentials [Enilij 
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Fig. 27. Ratios of (n,7) rates at T = 1 GK on stable and neutron-rich nuclei, obtained by using 
the HFB- 14 and d roplet model nuclear density distributions in the calculation of the optical 
potentials.CHSIM] 



where the disagreement between different models becomes larger. The main reason, 
however, is the lower relevant temperature because the differences in the rates due 
to the density distributions decrease with higher T, as the average transmission 
coefficients become less sensitive to the nuclear surface region. 

A further example of the impact of nuclear density distributions on calculated 
reaction S'-factors is shown in Figs. [151 [16] in the section on optical potentials above. 
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5.4.2.8. Deformation: Deformations are implicitly present when taking excited 
states from experiments or theory, or nuclear density distributions from theoretical 
models. There the problem lies in the fact that most results of microscopic models 
available for large-scale calculations assume sphericity. 

The HFM using transmission coefficients as described in ([64] and ((65|) cannot de- 
scribe reactions on deformed nuclei because it assumes ^ to be a good quantur n num- 
ber. The coupled-channel model has to be invoked for a rigorous treatment .^^^^^ 
It is computationally very expensive and thus not suited for large-scale calculations. 
Furthermore, the effective cross sections (jl9p require to include even more transi- 
tions than in standard, laboratory nuclear reactions. Fortunately, it has been shown 
that experimental data can be well described in a spherical HFM using an effective 
optical potential which is obtained by averaging over nuclear orientation.^^ This 
leads to a spherical potential with larger diffuseness. To compute this modified po- 
tential, the explicit inclusion of the nuclear deformation is required and is usually 
taken from microscopic or macroscopic-microscopic models. 

A deformation parameter may also enter the description of the NLD employed 
in the HFM calculation. But certain NLD descriptions include the deformation im- 
plicitly, like the one of Ref. |4] where the deformation is contained in the microscopic 
correction. 

The splitting of the GDR in deformed nuclei can also be accounted for phe- 
nomenologically by a dependence on a deformation parameter. This has the largest 
impact on the rates among the possibilities of the inclusion of deformation discussed 
here. 

Finally, fission transmission coefficients (see below) are also very sensitive to 
deformation parameters. Usually, the deformation is alrea dy included in an effective 
fission barrier, leading to double- humped fission barriers.!^ 

5.4.2.9. Width fluctuation corrections: The width fluctuation coefficients 
(WFC) defined in (|6ip impact the reaction cross section only closely around chan- 
nel openings, with a few keV to tens of keV. Contrary to isospin competition cusps 
(see above), the modify the cross section above and below the channel threshold.^^ 
Generally, they enhance the elastic channel and reduce the other open channels ac- 
cordingly to obey fiux conservation. Above the neutron threshold the behavior may 
be not so obvious because of the dominance of the neutron channel with respect 
to other channels. The transfer of strength from the dominant neutron channel to 
the elastic channel results in a marked reduction in competition with other exit 
channels and actually increases the strengths in those other channels. 

There are different ways to implement the calculation of the WFC, depending on 
the assumptions taken.!^ An explicit form for the W ^ requiring the knowledge of the 
level width distribution in the co mpound nucl eus, is obtained by assuming that all 
transmission coefficients are small .l^^^^^^^l^^ This approach has been widely used, 
even when the transmission coefficients were not small.^^ An alternative approach 
is to recognize that the main effect of the correlations is in the elastic channel. 
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with smaller effects on the other channels. This leads to a modification of the 
transmission coefficients and an additional factor applied in the elastic channel only. 
This is the more general HRTW method which has been used in most astrophysical 
applications.Q^^^ 200.201,97, Another implementation of the HRTW method gives 
a general formula wi thou t restricting assumptions for the additional factor but is 
complicated to apply.^^^ 

Overall, in astrophysical rates the impact of the WFC is small because of the 
low energies encountered. Even when the relevant energy window covers the neu- 
tron threshold, the difference when choosing one or the other description is barely 
noticeable because of the energy averaging taking place in the rate integration. The 
WFC are important, however, when trying to compare theoretical and experimental 
cross sections close to channel openings. 

5.4.2.10. Fission: Neutron-induced (and /3~-delayed) fission of extremely 
neutron-rich nuclides is important in determining the endpoint of the r- 
process aiid the a^^ nuclei produced by the fission 

process.l^^ '^ l '^^ l ^"^ ! ^"'^ ! ^"^ ! ^"^ ! ^"'' ! It can be included in the HFM by using an 
additional exit channel in (j64p describing the fission process. The fission trans- 
mission coefficient is calculated from the penetration probability through a fission 
barrier. Since most of the astrophysical fission occurs at energies below or close 
above the barrier, the resulting rates are very sensitive to the height and width of 
the fission barrier. Barrier predictions from various models show large differences 
and thus there are considerable uncertainties (reaching several orders of magnitude) 
in the resulting fission rates. These uncertainties have been explored in Ref. 12061 
and recommendations for comparative rate sets were given. The prediction of fission 
barriers remains a challenge to current microscopic models. 

Another issue concerns the fission fragment distribution. Earlier stu dies of the 
r-process have used simple fission barriers and assumed symmetric fission.l^SEZlSQS] 
Recent years have seen the advent of improved predictions using more so- 
phisticated statistical models and their results are being included in rate 
calculations! ^^^ l ^^^ l ^^^ l ^^^ l ^^^ l As for the barriers, considerable uncertainties may 
exist for extremely neutron-rich nuclei, though. 

On one hand, fission determines how far the r-process can synthesize elements 
and whether it could reach the region of long-lived superheavy nuclei. On the other 
hand, the fission rate together with the fragment distribution impacts the abun- 
dances of intermediate and heavy r-process nuclei. Intermediate r-process nuclei 
(including the rare earth peak) have a contribution from fission fragments.^^ZEOH 
If the fission timescale is short compared to the r-process timescale, fission cycling 
can occur whereby the fission fragments capture neutrons and follow an r-process 
evolution until they fission again. This exponentially enhances the final r- 

abundances. The number of possible fission cycles depends on the fission rates and 
thus is very sensitive to the fission barriers. The final abundance level is less sensi- 
tive to the fragment distribution but the distribution will determine the details of 
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which nuclei receive contributions from fission. 
5.4.3. A remark on HFM codes 

In principle, one would think that any code implementing the HFM should give the 
same result. This is obviously true when using the same implementation of the HFM 
and the same descriptions of the required input. However, both may differ among 
different codes. When quoting results it is therefore essential to not only always 
specify the exact version of the used program but also what selections regarding 
the properties described above have been made. 

There is a variety of codes which have been and are used in data evaluation. 
These usually focus on higher energies than astrophysically relevant. They also may 
include further reaction mechanisms which may not be relevant for astrophysical 
application (see Sec. I5.3.3|) . Finally, they use experimental knowledge (either di- 
rectly or by renormalizing theoretical results) or local parameterizations of nuclear 
properties. This way, high accuracy may be achieved locally for one or a few nuclei 
but no global prediction, essential for astrophysics, can be made. 

Astrophysical codes are especially written for global predictions of reaction rates 
and thus focus on low-energy cross sections. This includes using different internal 
numerics (see also Figs. [T71-t^for issues regarding the calculation of subCoulomb 
S-factors) but also different choices of the used input. Global treatments by global 
parameterizations or (semi-) microscopic models are preferred. Experimental infor- 
mation may be used where available to test those global approaches and, of course, 
to locally improve astrophysical reaction rates. However, a fair comparison of the 
global predictions of different astrophysical and other codes is only possible when a 
similar philosophy is used in determining their input (and by using the same exper- 
imental input where unavoidable). Most importantly, however, astrophysical codes 
directly account for the additional transitions required for the calculation of the 
true astrophysical reaction rate with thermally excited nuclei in the stellar plasma. 
They properly include the effective cross section defined in (fT9|) and thus implicitly 
use the correct weighting factors of transitions from excited states as derived in 
(gal). 

Early large-scale reactio n rate and cross section p redictions made use of a c ode 
developed at Caltech.HSESI xhe original codtPHEIll .^^^^g developed furthei^EI 
no further tables of reaction rates were published, ready for use in astrophysical 
reaction networks (see also Ref. |97] for a comparison of these early codes) . The 
work of Refs. [131 US] was not only important for astrophysical modeling but also 
in nuclear physics. Prior to these calculations, all hitherto experimentally studied 
reactions had featureless excitation functions and tests of the statistical theory were 
hindered by this. The large-scale calculations allowed to identify the cases suited 
to study competition between different channels.^^ ^^^^ 

Another influential development is the one of the SMOKER code.EZEH] 
went beyond the previous codes by including a more sophisticated calculation of 
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the transmission coefficients in all channels by explicitly solving the Schrodinger 
equation with optical potentials and including several new global parameterizations 
of the further inputs. A set of neutron-capture rates for r-process nucleosynthesis 
was published in Ref. |26] and complete sets of neutron-induced as well as charged- 
particle induced reactions calculated with this code became included in the first 
version of REACLIB, a library of theoretical and experin iental cross sec tions which 
can directly be used in astrophysical reaction networks P^ I "^'*^^ ! "^'*^'^ ! "^'*'^ ! The REA- 
CLIB format, using fits of reaction rates to a function of 7 temperature-dependent 
terms (see, e.g., Re f. IHl), has become a standard widely used in the astrophysics 
communityl ^ l ^^ l ^-^ ^ 

A series of codes - some closely related, other only loosely based on SMOKER 
- has appeared since then. The code NON-SMOKER extended the functionality 
by not only updating the inp ut data b ut also including an improved, global NLD 
description and isospin effects. I^^HHEHI Although it already included many possibil- 
ities for descriptions of nuclear properties to use, often the term "NON-SMOKER 
calculations" is used synonymously for the extended tables in Refs. [9l |34l [220] cal- 
culated with a chosen input set. In a parallel development, the code MOST also 
updated t he inp ut physics and provided a different selection of treatments of nuclear 
properties! -*^^^ ! Tables of reaction rates (but no fi ts in REACLIB format) were pro- 
vided online for several versions of the code.^^^ The newly w ritten , but still closely 
related, NON-SMOKER^^^ code included several changes.!^ Apart from sev- 
eral updates of included nuclear data, also the internal numerical calculations were 
modified, the isospin suppression was improved, additional choices for microscopic 
and macroscopic predictions of nuclear properties were offered. The innovative web 
interface allows access to the code from anywhere through a simple web browser. 
Additional switches can be set and different nuclear properties provided in an op- 
tional input file which is uploaded to the server running the code. The resulting 
cross sections and reaction rates are immediately displayed. Certain nuclear prop- 
erties, e.g. optical potentials, cannot only be uploaded as data but also as formulae 
because the code includes a simple equation parser. NON-SMOKER^^^ has been 
used for the astrophysical analysis of a large number of experimental results, espe- 
cially for p-process nucleosynthesis (see, e.g., Refs. I141[ I142[ [22l [T43l 11441 11451 and 
references therein). Its development has been frozen at version v5.8.1w but it is still 
available and used for calculations. 

Most recently, the code SMARAGD (Statistical Model for Astrophysical Reac- 
tions And Global Direct reactions) continues and extends the development initiated 
with NON-SMOKER ^^^ of a u ser- friendly, easily extensible code tailored for astro- 
physical reaction rates .'^SlMlMli'he code is written completely in FORTRAN90/95 
(with exception of the routines handling the web interface and the function parser, 
which are written in C) and has a modular structure, making changes easy. The 
latest nuclear data can also swiftly be included through files or web downloads. 
Internally, the numerics and solvers for the Schrodinger equation have been im- 
proved to be more accurate at low, subCoulomb energies and to consistently also 
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calculate direct processes as required for astrophysics (see Sec. I5.5.ip . Recent de- 
velopments in parameterized or microscopic predictions of global nuclear properties 
(masses, NLDs, density distributions, optical potentials, photon strength functions) 
have been included. It also uses the modified HFM discussed in Sec. 15.3.31 Reaction 
rates are provided in tabular form as well as in the REACLIB format through an 
automated fit routine. Future versions will include multiple particle emission, fol- 
low 7-cascades explicitly, and allow the calculation of fission rates. Code versions 
below vl.Os are not public, those below v2.0s do not include direct reactions, yet. 
Direct reactions are included in a number of ways, as discussed in Sec. l5.5l below. At 
later stages of the code development it is planned that users may upload modules 
providing nuclear data, numerical methods to compute required properties, or even 
altering the functionality of the program. However, before these improvements in 
versatility are made accessible, a new large-scale calculation will provide a new set 
of published reaction rates between the driplines, intended to improve on and su- 
persede the NON-SMOKER rate s'^ ' ^ which currently are used by the majority 
of astrophysical modelers worldwide.l^ 



5.5. Direct reactions 
5.5.1. General remarks 

In its most general definition, the term "direct reaction" includes all processes di- 
rectly connecting the initial and final states of a nuclear reaction without formation 
of an intermediate compound system. This includes elastic scattering as described 
in the optica l model, and inelastic scattering which predominantly excites collective 
states.f^Sl^ The latter includes Coulomb excitation which has been found to be im- 
portant in heavy ion collisions due to the high Coulomb barriers involved. ^^^IZZlllS] 
In astrophysically relevant reactions, especially with a-particles, energies may also 
be close to or below the Coulomb barrier and Coulomb excit ation may also become 
important, depending on the structure of the target nucleus.^^SH 

Here, we focus on direct reactions when some (if it is a stripping reaction) or 
all (if it is a capture or charge-exchange reaction) nucleons of the projectile are in- 
corporated in the target nucleus. In a pick-up reaction, one or more nucleons from 
the target nucleus are added to the projectile to form the ejectile, again in a direct 
manner. Pick-up and stripping reactions are subsumed under the term "transfer 
reactions" . In contrast to the HFM model, direct reactions excite only few degrees 
of freedom because most of the nucleons included in the system of target nucleus 
plus projectile remain spectators. A nucleon of the projectile reaches its final state 
without sharing any energy with any of the other nucleons present and the ex- 
cess energy is emitted as a discrete photon carrying the energy difference between 
initial and final state. Direct reactions can be identified experimentally because 
of their angular dependence of the differential cross sections, being peaked in for- 
ward direction. Direct processes are also faster by at least 5 orders of magnitude 
than compound reactions, with reaction timescales of the order of 10^^^ s. This is 
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Fig. 28. Relation between di rect n eutron capture and compound capture in the HFM as function 
of neutron separation energy.GSS The result for Ti, Cr, Fe are taken from Ref. 12351 the ones for 
Sn from Ref. 12361 and the ones for Pb from Ref. 1651 



comparable to the time the projectile requires to cover a distance of the size of a 
nucleus. Therefore, direct reactions are important at high projectile energies when 
compound formation is disfavored. 

Although the notion of direct processes was inspired originally by angular distri- 
butions of low-energy reactions, it was assu med fo r a long time that higher energies 
are the domain where they are dominating .l^^l^ In resonant reactions at lower en- 
ergy, it is sometimes necessary to include a non-resonant background (which may 
show interference with resonances) but experimentally it is often difficult to dis- 
tinguish between a direct component and contributions from tails of resonances. 
However, in systems with low NLD, and thus widely spaced resonances, direct re- 
actions become important even at astrophy sically low intera ction energi es because 
compound formation is suppressed.^^^ 226 227 | 228 | 2 29 230 23U2 32 | 233 | r^j^jg ^^^^ 

applies to intermediate and heavy nuclei far off st ability, e.g., for neutron capture 
in nuclei with low neutron separation energy.l^llSSI xhe direct capture cross section 
can become considerably larger than the compound cross section. Figure [28] com- 
pares (n,7) cross sections at 30 keV from direct capture and from HFM for a number 
of isotopes. With decreasing neutron separation energy, the direct component plays 
an increasingly important role because the compound nucleus is formed at lower 
excitation energy and thus also at lower NLD. Elements with inherently low NLD, 
such as Sn, show a larger direct contribution for all isotopes. Also for nuclei at shell 
closures the NLD is low and the importance of direct reactions enhanced relative 
to the HFM. ■ Similar considerations may apply to proton-induced reactions on 
proton-rich nuclei. 

Direct processes are not only important to be included in reaction rate predic- 
tions. All possible processes (elastic and other direct ones, compound-elastic and 
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compound nucleus reactions) have to be taken into account in the analysis of scat- 
tering and reaction data when extracting optical potentials. The elastic scattering 
cross section cTg (which is a direct process), the one for reactions (Tr, and the total 
(Tt are related by 



where we distinguished between elastic scattering at the optical potential cr°P* and 
compound-elastic scattering (Tg°™P. The cross section CTx includes all inelastic pro- 
cesses, i.e. reactions. Often, one reaction mechanism is dominating by far and then 
(Tx can be identified with the cross section for that mechanism, e.g. the compound 
cross section cr™™P for resonant processes (including the HFM) as discussed in 
Sees. 15.21 and 15.3.21 But this is not always the case as - depending on the nucleus 
and the projectile energy - there may be additional mechanisms contributing to Ux 
in some cases. Therefore, using an optical potential derived from scattering implies 
that the absorption term is due to some reaction(s) but does not define the reaction 
mechanism(s). Using such a potential in a pure HFM implicitly assumes that the 
missing flux from the elastic channel is due to the compound mechanism only. This 
may not be appropriate when direct processes are non-negligible (this comprises 
direct reactions as discussed below but also Coulomb excitation at low energy) and 
will require a modification of the optical potential depending on which mechanism 
is to be studied. This is also implicitly contained in the idea of the modified HFM 
briefly discussed in Sec. 15.3.31 Earlier calculations of direct neutron capture have 
made use of a simple hard sphere capture model (see Sec. l5.5.3|) in a combination of 
direct and HFM capture, not just to simplify calculations but also because it also 
allows the assumption that the contribution to direct capture potential absorption 
due to the tail of distant resonances is already included in the statistical model 
averaging.!^ 

Astrophysical rates can be calculated from cross sections by applying ([5]), ((HI), 
and (|25|) . regardless of the reaction mechanism. Since each discrete transition ap- 
pearing in a direct process obeys the reciprocity relation (|17p . a similar effective 
cross section (|19|) can be derived as for the compound case, resulting in the same 
weighting factors (|42|) of transitions from the excited states. Also the same reci- 
procity relations (|27p and ([28]) apply, provided that thermal population of the states 
in all participating nuclei is valid. 

In the following only a few methods are outlined which have been used to cal- 
culate direct reactions for astrophysics. Only the basic equations for the reaction 
from one initial state to one final state are given but the actual rate equations can 
straighforwardly be obtained by using the methods described in Sees. 14.11 and 14.31 
It should be noted that here not only light targets are implied but also interme- 
diate and heavy ones (see Fig. [28]) . for which microscopic models are not feasible 
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(see Sec. l5.1|) . Furthermore, low projectile energies are implied as these are required 
for calculating astrophysical reaction rates. Nevertheless, direct reactions at higher 
energies can be used to extract certain properties, such as spin assignments and 
spectroscopic factors (see, e.g., Ref. I237p . of stable and unstable nuclei which are 
required for the calculation of the cross sections and the rates. 



5.5.2. DWBA 

Direct transfer reactions can be treated by solving the time-independent 
Schrodinger equation with optical potentials in the entrance and exit channels. A 
simple implementation of this is the Distorted Wave Born Approximation (DWBA) . 
The differential cross section for the one-nucleon or cluster transfer a + A ^ b + B 
w ith a — x = b, A + x ~ B for Aa < 4 and = 1 or 3 is given in zero-range DWBA 

daf'" _ uiAamBbkBbgB s;^ at '^^^3^^^ l-7^\ 

with the zero-range normalization constant Nq. The reduced cross section without 
spin-orbit coupling is given by 

a,sj{0)^Y.K\' , (76) 

m 

with the reduced transition ampHtude 

til, = ^ / * (^kst, ue.ir) [^'Yrir^] * x^^J (kAa^r) dr . (77) 

As before, the quantities mAa, fnsb and kAa, ksb are the reduced masses and 
wave numbers in the entrance and exit channel, respectively. The orbital angular 
momentum quantum number i, the spin quantum number s, and the total angular 
momentum quantum number j refer to the nucleon or cluster x bound in the residual 
nucleus B. The spectroscopic factor and the isospin Clebsch-Gordan coefficient 
for the partition B = A + x are given by C and Sij, respectively. The optical 
wavefunctions in the entrance and exit channels are given by x^~^^ and the time- 
reversed solution x^""* ■ The bound state wave function is denoted by Uf j and the Y/" 
are the usual spherical harmonics. Expressions similar to the above are obtained 
when the finite range of the interaction potential is taken into account.^^ 

Important for the successful application is to keep the number of open pa- 
rameters as small as possible. For this reason, folding potentials (see the para- 
graphs on optical potentials in Sec. I5.4.2|) were used in many astrophysical appli- 
cations of the model, with A either determined from scattering data or from global 
dependences. This leaves the spectroscopic factor S which is usually de- 

termined by simply comparing the calculated magnitude of the differential cross 
section to measurements. This is mainly do ne w ith (d,p) or (d,n) reactions at ener- 
gies above the astrophysically relevant ones.^^^ The partial width F^ appearing in 



June 10, 2011 0:11 WSPC/INSTRUCTION FILE rauscher 



78 T. Rauscher 



(|48p can be related to spectroscopic factors for a particle z in a state fj. 





offers a different experimental access to spectrosco pic facto rs (see Sec. 15.5.51 for 
further methods to determine spectroscopic factors). ^^SSHMl 

Further required input includes, of course, masses or separation energies and 
nuclear spectroscopic information which have a similar importance as for the HFM, 
discussed in Sec. 15.4.21 

In the absence of experimental data, the spectroscopic factor can be calculated 
microscopically from the overlap between initial and final state wave functions, e.g., 
in the shell model. '''^ ''^ ^'^'^ '^•^ ■^^'^ ^'^ However, there is some amb iguity because this 
overlap is not well defined in different microscopic approaches .^^^^^ It has been 
shown that the spectroscopic factors for depositing or picking up a single nucleon 
are related to t he occup ation factors of the participating quasi-particle states (see 



As in the case of the HFM, many codes have been used for DWBA calculations 
over the years, especially for the analysis of data at intermediate and high energies. 
The code TETRA has b een written especially for application at astrophysically 
relevant low energies .l243jj^ j^g^g |-)ggjj applied successfully to astrophysically relevant 
reactions with light and intermediate target nuclei (see, e.g., Refs. 12251 12261 



Is it necessary to go beyond the DWBA? There are three fundamental assump- 
tions contained in the DWBA treatment!^ 

(1) The reaction proceeds directly from initial to final state and all particles except 
the transferred one(s) remain unaffected spectators. 

(2) The wave function for the relative motion between the reactands is assumed to 
be correctly described by the optical potential. 

(3) The reaction is assumed to be sufficiently weak to be treated in lowest order. 

To relax the first two assumptions, the coupled-channel Born approximation 
(CCBA) was introduced.^^SESl Under rare circumstances the transfer amplitudes 
may be large and the third assumption h as to be relaxed. This leads to a full 
coupled-channels treatment for the reaction.l^^^^SEH 

One has to be aware of the fact that the relevant energies remain low for as- 
trophysical reaction rates, also because transitions from excited states contribute 
considerably. This is contrary to what one is used to in the investigation of reactions 
proceeding at several tens of MeV. Due to the low energies involved, the reaction 
channel is weak (compared to, e.g., elastic scattering) and the third assumption is 
valid. The usual concern with the second assumption is that the optical potential 
also has to describe well the wave function even deep in the nuclear interior. The 
deep region, however, is crucial for reactions at higher energy whereas at the low 
astrophysical energies most contributions to the overlap integrals stem from regions 




mm USD mi IMl HH IMZI and references therein). 
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close to the surface of the nucleus or even from outside of the nuclear radius. As 
long as these regions are described well by the optical potentials, the DWBA should 
work. The first assumption implies that either no indirect processes exist or that 
they can be treated separately (incoherently) as was suggested by the above, sep- 
arate discussion of compound reactions and other mechanisms. In the HFM it is 
assumed that interference terms cancel and thus also interference with direct reac- 
tions should cancel on average. Interference with isolated resonances can be treated 
explicitly by adding an interference term, e.g., between the S-factor of the direct 
reaction ^direct and the one of a Breit-Wigner resonance Sbw (see Sec. 15.21) 

S = ^direct + Sbw — 2(v^5^iircct5W') cos Winter , (79) 

where 

Winter = arctan f ^ , ° ' j (80) 
is the energy-dependent, relative phase shift.^^lS] 



5.5.3. Direct capture 

A potential model can also be used to calculate direct capture (DC). Although 
microscopic models are an alternative for light systems (see Sec. l5.ip . a DC potential 
model has the advantage that it can be applied also to heavier nuclei. The DC cross 
section for a particular transition is determined by the overlap of the scattering wave 
function in the entrance channel, the bound-state wave function in the exit channel, 
and the electromagnetic multipole transi tion oper ator. 
The DC cross section is then given bjISSSESQ] 

= [dn2(^^^)(^y^ |tM.M.M.,.r . (81) 

J V J \kAaJ 9Aga 

The polarisation a of the electromagnetic radiation can be ±1. The wave number 
in the entrance channel and for the emitted radiation is given by kAa and k^, 
respectively. 

The multipole expansion of the transition matrices TMAMaMB,cr including elec- 
tric dipole (El) and quadrupole (E2) transitions as well as magnetic dipole (Ml) 
transitions is given by 

tMAM^MB,a- = tfl^M^MB,a dsai^^) + tA'lAMaMB ,cr '^L (^) + ^MiM„MB ,(t 4(t(^) ■ (82) 

The rotation matrices depend on the angle between kAa and k-y which is denoted 
by e, where S ^ Ma + Ma - Mb- 
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Defiiiing 



C(E2) 



Aa Aa j 



Ma 



(83) 



we can write for the transition matrices for the electric dipole (E£ = El) or 
quadrupole (E£ = E2) transition 



= i'" (^a Sa Ma \ ja Ma){]b Mb - Ma I A Ma I Ib Mb) 



laja 

X (£ 5 jb Mb - Ma \ ja Ma) C{EC) la h h 

X {hOCO\laO)W{ChjaSa;lajb)lF,il,;l^,^ ■ (84) 

In the above expressions the quantum numbers for the channel spin in the entrance 
channel and for the transferred angular momentum are denoted by ja and /q, respec- 
tively. The quantities I a, Ib and Sa (Ma, Mb and Ma) are the spins (magnetic 
quantum numbers) of the target nucleus A, residual nucleus B and projectile a, 
respectively. 

For magnetic dipole transitions (M£ = Ml) we obtain 



t 



MC 

MAMaMB,(r 



J2i''<j\ {la Sa Ma \ Ja Ma){jb Mb - Ma I A Ma I Ib Mb) 



X {IS jb Mb -Ma \ jaMa) 



( Za 










mi J 



+ 2^1a{-^)''-'- Sa Jb VSa{Sa + l) G{1 Sa ja U Sa Jb) 
- {la Sa Ma \ Ja Ma){ja Ma Ia MB-Ma \ Ib Mb) 

X {Ia Mb- Ma I 5 \IaMa) 

] ( he ^ 

X MA 5j^j^ V {Ia + 1)/^A 



2mpC^ 



r / rMl 



(85) 



where G is the Racah coefficient, the are the magnetic moments and rrip is the 
mass of the proton. 

The overlap integrals in (|84| and (j85]) are given as 



/ druNLj{r)0'^^{r)xu,^{r) 



EC, 



(86) 



for the electric dipole (EC — El) or quadrupole (EC = E2) transition, and by 

druNLj{r)0^\r)xi^,^{r) (87) 



lbjblB;laja 
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for the magnetic dipole transition (M£ — Ml). 

The radial part of the bound state wave function in the exit channel and the 
scattering wave function in the entrance channel is given by UNLj{r) and Yi (r ), 
respectively. The radial parts of the electromagnetic multipole operators are 

'"(r) = — - [sin p - 



-iMl/ 



OEi(r) 



2/5 
3 

15 



+- p cos p\ , 
2) sin p + 2/5 cos p\ r 
- 12)sinp+ (12 



p'^)pcosp\ 



(88) 



In the long wave-length approximation 
quantities reduce to 



applicable as long as p ■ 



fc-yT <c 1 - these 



O^i(r) ~ 1 



(89) 



Usually, only the dominant El transitions have to be taken into account. Possible 
exceptions are captures far from stability with very low reaction Q-values because 
for these cases no final states may be energetically accessible through El transitions. 
However, because the astrophysical reaction rate involves summing over transitions 
originating from excited states, a larger spin range may be available and El (from 
excited t arget states) may again dominate. For El transitions, the above expressions 
reduce td^^ 



'El 



IQtt 

E 



iff 




kAafic J \n/ gag A 
{2Jfs + 1) (2Ja + l)max(£a,^^) 



1 C 



(90) 



The coefficients aj are calculated in LS coupling to 
a^j = gA{2I + 1) {2Lb + 1) {2Sb + I) 



I La Sb\ ( I La Sb\ 
Lb Ib ip},\lA 4 /a Je ■ ^ ' 

In the above expressions, the energy of the emitted photon is Ej. The orbital and 
total angular momentum quantum numbers of the nuclei in the entrance and exit 
channels are ia, Ja, and Jg, respectively. The spin quantum number, the orbital 
and total angular momentum quantum numbers are characterized by 5, L and /, 
respectively, with indices a, A and B corresponding to the projectile, target and 
residual nucleus, respectively. The notation {• • • je stands for the 6j symbol. The 
radial wave functions in the entrance and exit channels are given by Xa and Ujs, 
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respectively. The spectroscopic factor and the isospin Clebsch-Gordan coefficient 
for the partition B = A + a are given by C and Si^ , respectively. 

The DC potential model has been successfully applied to many reactions 
with light targ et nuclei at astrophysical energies, e.g., for ''Be(p,7)^B and 
^Li(n,7)8Li.l2Sl For further calculations see, e.g., Refs. |M1 IMl IMl 1221 [ISl 
[2481 [2^ [2541 [2551 [2561 [2571 [2581 [2591 and references therein. Many calcula- 
tions were performed with the DC potential model code TEDCA, which is 
tailored to treat low-energy reactions of astrophysical interest .^26Q| Again, fold- 
ing potentials (see Sec. I5.4.2p were the key to reduce the number of parame- 
ters. Spectroscopic factors were taken from experiment or shell model calcula- 
tions. It was also used t o extend calcu lations to intermediate and heavy target 
nuclei.tl^O 229 233 235 236 24()I247I26T12 62 263 264 265 266 267 268 

Other astrophysical calculations have not made use of the full DC equations 
but used simplifying assumptions. Resonant and DC rates on the proton-rich side 
based on hard sphere scattering wavefunctions in the entrance channel - with res- 
onance properties, final states, and spectroscopic factors taken from shell model 
calculations - were provided for target nuclei with mass number 44 < A < 63.1^ In 
these calculations, however, the resonant part exceeds the direct part of the cross 
section by several orders of magnitude. Astrophysical neutron capture on stable 
and neutron-rich nuclei was calculated in the hard-sphere model for El capture by 
Refs. [Ml UMl W!^JFn[J Fn\ In this model, the El neutron capture cross section 
can be written a3^2^^ 



(92) 

where rhard is the hard sphere radius and the multiplicity ^ is the number of incident 
channel spins which can lead to the same final state with spin J"^. It is ^ = 1 for 
J'X = 0, or for J'X^Q and J" = J'X± 3/2. The value £, = 2 applies for 7^ 
and = ± 1/2. The quantity l'^ is the orbital angular momentum of the final 
bound state v and Sj^ is the spectroscopic factor of this state. The dimensionless 
parameter Y is given by 



V> ^hard 

h 

The advantage of this approach is that no explicit wave functions for scattering and 
bound states are required. On the other hand, the correct overlap of the wave func- 
tions may yield more accurate cross sections, especially for low projectile energies 
when considerable contributions to ()86|) are coming from far outside the nuclear 
radius. Ref. [61] showed the importance of direct neutron capture in the r-process 
based on the hard-sphere model. 

For the validity of the low-order potential model approach for DC while neglect- 
ing higher-order processes, similar arguments can be made as were presented for 
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Fig. 29. Comparison of direct neutron capture cross sections of Sn isotopes calculated with input 
taken from a Hartree-Fock-Bogolyubov model (HFB, Refs. [2741 1275|l . a Relativistic Mean Field 
model (RMFT, Refs. 12761 12771 1. and a semi-microscopic model using folded Yukawa potentials 
(FY, Refs. [2781 [2791 [280)l . For details, see Ref.[65l 

the DWBA towards the end of Sec. \5T2\ 

5.5.4. Sensitivity of DWBA and DC in astrophysical calculations 

For astrophysics, a large number of rates for highly unstable nuclei have to pre- 
dicted. Similar to the uncertainties of the HFM discussed in Sec. 15.4.2] the DWBA 
and DC predictions are sensitive to certain inputs, such as masses (Q- values), ex- 
cited and bound state properties, spectroscopic factors, optical and bound state 
potentials. When folding potentials are used, a (weak) dependence on nuclear mat- 
ter density distributions appears additionally. The discussion of these quantities 
in Sec. 15.4.21 also applies here. However, some of the quantities appearing in the 
HFM are sums over individual transitions and averaged quantities whereas indi- 
vidual transitions are determining the direct reaction cross sections. Thus, direct 
reactions are more sensitive to nuclear properties impacting those individual tran- 
sitions, including spins, parity, energy of bound and excited states, Q-values, and 
spectroscopic factors. Due to the angular momentum barrier and the low projectile 
energies, low partial waves contribute most to the El cross section, i.e. s-waves 
when initial and final states have differing parity, p- waves when they have the same 
parity. Assuming equal spins and parities, low-lying states are contributing more 
than higher ones. The strong sensitivity to individual transitions, however, may be 
reduced in the astrophysical rate, when using the weighted sums ([8|), (|25|) . (|42|) over 
transitions from excited target states to available final states. 

For neutron-rich Sn and Pb isotopes, predictions for 30 keV neutron capture 
employing input from different microscopic or semi-microscopic approaches were 
compared in detail.l^The DC calculation was tested in comparison to experimental 
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Fig. 30. Same as Fig. 1291 but for Pb isotopes. 



data for ^°^Pb(n,7)^''^Pb. In this case, a discrepancy between the data obtained in 
an activation measurement and the data from a high-resolution resonance counting 
experiment was resolved by showing that the difference is due to the DC contri- 
bution to the cross section which is only included in the cross section from the 
activation measurement. For very neutron-rich isotopes of Sn and Pb it was found 
that the resulting cross sections differ by orders of magnitude with the different 
inputs (see Figs. [29l [30)) . This is mainly due to the sensitivity of the cross section 
to the predicted location of the low-spin bound states with respect to the neutron 
separation energy.iSS The DC calculation can be tested for ^'^^Sn(n,7 ) because there 
is experimental information on the bound states in i33gj]^[237l281] This nucleus is 
predicted to be close to or directly in the r-process path.^^ An independent cal- 
culation confirmed the original work.l^Sl Fortunately for astrophysics, very similar 
cross sections (within a factor of 3) are computed for this reaction with input from 
the different microscopic approaches. Unfortunately for nuclear physics, the reac- 
tion is not a good case to select a preferred microscopic model for the same reason. 
Spectroscopic data for neutron-richer isotopes would be necessary. Furthermore, 
spectroscopic factors were set to Unity in these calculations. This is a good as- 
sumption for the states in ^^^Sn (and was recently confirmed by Ref. I237P but is 
not valid in mid-shell. However, the uncertainty introduced by the predicted bound 
state energies exceeds by far the one introduced by the spectroscopic factors when 
only particle states are considered and hole states are neglected. '^'^ The latter have 
very small spectroscopic factors and are negligible because they involve a reordering 
process in the final nucleus (see also Ref. 12351 for another example of a reordering 
process in ^^S and its dependence on deformation). 
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5.5.5. Averaged direct capture 

Regular DC cross sections are obtained by summing over all allowed transitions to 
energetically accessible final states, 



'DC 



(94) 



where each summand contains the appropriate spectroscopic factor. Following the 
derivations in Sec. 14. 1[ the astrophysical rate will contain an effective cross section 
and a partition function. Similarly to the HFM, it was suggested (Refs. 12701 12711 
I272[|282)) that the sum over final states may be (partially) replaced by an integration 
over the NLD in the final nucleus pf, 



'DC 



(E) 



H^DC 



+ 



+ / 2^ pf{Ef,jf,TTf)a^-^^{E,Ef,jf,TTf,SEfjf.f)dEf . (95) 

If jfTTf 

Similarly, the summation over initial states p, in (|19p can be (partially or fully) 
replaced by an integration over the NLD in the target nucleus pi, 



O^DC 



'DC 



^5proj(T) 



E 



p,{E\J\7:')a 



DC 



{E, E\ J\tt\ SE'jiT^'EfjfTTf] 



dE' 



(96) 



Through the NLD, the cross sections and not only include transitions to 
discrete final states but also "average" transitions to states described by the NLD. 
Therefore I call this averaged direct capture (ADC). As in the HFM, the ADC cross 
section may only include the ground state transitions and NLDs above the ground 
state when all excited state properties are unknown. The advantage of this is that 
the sensitivity to the location of discrete states relative to the projectile separation 
energy (as seen in Ref. I65p is washed out and the change in cross section from one 
isotope to the next is smoother. Of course, this may not properly describe nuclear 
structure details in the cross sections when only few states are available but it 
can give a more appropriate estimate of the magnitude of the rate for astrophysics 
(which involves an averaging over a relevant energy range and includes more tran- 
sitions than the laboratory rate) than relying on just a single (semi-)microscopic 
approach. 

Obviously, the ADC cross section will be sensitive to the NLD. Contrary to the 
HFM, where the NLD mostly impacts the 7-widths and has its strongest effect at 
several MeV excitation energy, captures to low-lying states (assuming spins and 
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parities are favorable) are dominating and thus the NLD at low excitation energy 
will be relevant. This is why it is important to use a proper spin- and parity- 
dependent NLD description. It was found recently, however, that thermal excitation 
of target nuclei reduces the sensitivity to the parity dependence in the NLD for 
astrophysical rates in the HFmISSI We can expect a similar effect for direct capture, 
although the number of possible transitions is more limited-^SQ] 

Appropriate spectroscopic factors are a further important ingredient in ADC 
calculations. In the integral of ([95)) these are averaged spectroscopic ia.cioYS Ssf jf-Kf , 
describing the average overlap between the initial state ji and the final bound 
states with given spin and parity at an excitation energy . The doubly averaged 
spectroscopic factors SEijiTr'Ef jfTrf appearing in the integrand of ([M)) are even 
more complicated, as they involve the average overlap between all initial states 
with given spin and parity at excitation energy and all final states. Currently, 
the only spectroscopic factors to be found in literature (from experiment or theory) 
are for transitions connecting the ground state /i = of the target nucleus with the 
final states. Spectroscopic factors for transitions from excited states have yet to be 
calculated. They are needed not just for the averaged DC model employing NLDs 
but also for the regular DC model when applied to compute astrophysical reaction 
rates (see also the discussion of the stellar enhancement and the effective weights 
in SedOJ. 

Spectroscopic factors for one-nucleon capture (or transfer) on a target in the 
ground state can easily be computed fro m the occupation nu v"^ as calculated, 

e.g., from BCS or Lipkin-Nogami pairing.l ^^ l '^'^^ l ^^ l "^^^ ! "^^^ ! Then the spectroscopic 
factor for putting a nuclcon in state j with spin Jj is just 

S, = l~ v] (97) 

for a target nucleus with an even number of nucleons of the same type as the 
projectile. The occupation probability is the one of the target nucleus. In a chain of 
linked reactions, the total processing efficiency is given by the slowest reaction(s). 
Therefore, when considering sequences of capture reactions, e.g., in the s-, r-, rp- 
processes, the rates on such target nuclei with even nucleon number will be the 
slowest (and their reverse photodisintegrations the fastest) and thus they will have 
the largest astrophysical impact. 

When the target nucleus has an odd number of nucleons of the projectile type, 
the spectroscopic factor is 

S, = (2J,^\)v] . (98) 
The occupation probability w| is always taken in the nucleus with even number of 



nucleons. In ([98|1 this is the final nucleus of the reaction. The expressions for the 
extraction of a nucleon from a given state follow from the fact that the time-reversed 
reactions have to show the same spectroscopic factor. The occupation probabilities 
can be calculated from m icroscopic theo ry, e.g., using BCS or Lipkin-Nogami pairing 
on a single-particle basis .1^219012351279] 
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Fig. 31. Comparison between experim enta l and averaged (Eg. llOll spectroscopic factors as func- 
tion of excitation energy for iSSxe+n.EISl 



In the absence of calculated spectroscopic factors several different approaches 
have been used in the past. Often, spectroscopic factors were set constant to 1.0 
(e.g., Refs. ins 12331) or to 0.1 (e.g., Refs. UHl IMS)- These values can already 
be seen as averaged spectroscopic factors iS^/j/^/ for low-lying particle-states. 

Useful for the application to the ADC approach is the construction of an exci- 
tation energy-dependent, averaged spectroscopic factor to be employed along with 
the NLD. As an example for this, in studying neutron capture on the astrophys- 
ically important nucleus ^""^Ti the DC component was estimated implementing a 
distributed spectroscopic strength.I^^ In this case, transitions to l/2~ and 3/2~ 
states in ^^Ti are dominating. A distribution (j)^ of the l/2~ and 3/2" strengths 
was assumed, reaching from the (experimentally known) location of the lowest l/2~ 
and 3/2~ state, respectively, to the neutron separation energy. With this "smear- 
ing out" of the states and of the strength and due to the E!^ dependence of the El 
transition probability, then the calculation is reduced to computing a transition to 
an effective bound state with full spectroscopic strength at an energy of 

Abound = , (99) 

where the average transition energy is given by 

Such an approach accounts for the uncertainties in spectroscopic strength and lo- 
cation of excited states and can be viewed as a zeroth approximation to the ADC. 

Another suggestion for the functional form of the energy-dependence of the aver- 
age spectroscopic factors was made in Refs. I270T 12711 [?7^ The spectroscopic factors 
describe the overlap between the antisymmetrized wave functions of target -|-nucleon 
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and the final state. The number of final state configurations increases with increas- 
ing excitation energy E'^ and the overlap of initial and final state wavefunctions 
decreases. Thus, also the spectroscopic factor decreases. In a simple approxima- 
tion, the energy dependence of the spectroscopic factor for single nucleon transfer 
can be parameterized by a Fermi function with 

S = h.rUE^) = \^ (101) 

1 -I- e A* 

and the parameters E* ,A* . This is motivated by the excitation-energy dependence 
of the occupation probabilities. Figure shows how well averaged spectroscopic 
factors with the functional dependence (I101|) compare to experimental ones for 
I36xe-Hn.ll711 

There is also a connection between single-particle spectroscopic factors Ssp and 
the partial resonance widths F = F^, F, . . . appearing in the BWF (see Sec. 15. 2p 

T = 2PeC^Sspel^^ , (102) 

where C is the isospin Clebsch-Gordan coefficient, m the reduced mass of 
the system nuclcus+particle, and 0q is the dimensionless single-particle reduced 
width.l^ •^'^^ ■^^ The penetrability Pi for the relative angular momentum i 
can be expressed in terms of the regular and irregular Coloumb wavefunctions Fi 
and Gi 

In literature, a value of 9q — 0.6 is often assumed for an average single-particle 
reduced width. Comparing ([78]) and (jl02|) . it can be seen that an average 6q can 
be calculated from solution s of the radial Schrodinger equation (giving F) with 
an optical potential.l^^^^^^ Since (F) = Ttrans/(27rp), there is a direct connection 
to the average transmission coefficients Ttrans — , T, . . . appearing in the HFM 
(see Sec. 15.3.2]) which can be used to consistently derive the average combined value 
(S) = (iSsp^o). The averaged spectroscopic factors for transitions from excited target 
states in can be estimated in this way. 

The integrals appearing in (^51) and (^51) always contain products p{S) of NLD 
and averaged spectroscopic factors. With the relation between strength function 
Sf — 7trans/(27r) and reduced width (see, e.g., Ref. I288P we can derive 

27rTOrnuc^ 

P{S) = — = ^trans (104) 

and thus (partially) eliminate the NLD within the integrals as it is implicitly con- 
tained in the or Ttj-ans- The strength function is only defined for resonance 
states above the projectile separation energy. For estimating transitions to bound 
states it was suggested to construct an "internal strength function" s^-n t , extending 
the regular strength function below the separation energy.EZQEZlIZISl ^jj^jg jg j^q_ 
tivated by the relation ()104|) and the fact that the spectroscopic factors describe 
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Fig. 32. Sums of the contributions to the total direct capture cross sections with the averaged DC 
model and the standard DC for ^'^^Xe(n,7) as function of excitation energy of the final nucleus.EZS 



the structure of bound states as well. The ADC was investigated in Ref. '272' for 
the reactions -'^^^Xe(n,7)-'^^^Xe and "'^*''^Sm(n,7)-'^*^Sm, using (|92l) with an internal 
strength function. The ADC cross section (for the target in the ground state) was 
compared to a calculation summing over transitions to bound states in the stan- 
dard potential model (as described in Sec. 15.5.3^ and with experimental data. The 
energy-dependence of Sint was chosen very similar to the one found for the average 
spectroscopic factors Sint = C'*/FGrini, with independent parameters C* , E* , A*. 
The parameters were determined by requiring Sint = s/ at the neutron separation 
energy. Simultaneously, it was required that the ADC cross section integrated up 
to the excitation energy of the last included state yields the same cross section 
as obtained with the standard potential model. Although this is not suited for a 
prediction, it can be used to assess the validity of the assumptions. Including only 
s-wave neutrons, the experimental cross sections were reproduced within 25% for 
"'^^^Xe(n,7) and 2% for ^*^Sm(n,7). Figure [32] displays a comparison between the re- 
sults from the averaged direct neutron capture and standard direct neutron capture 
for i36Xe. 

For predictions across the nuclear chart, the spectroscopic factors and/or in- 
ternal strength functions can be obtained from optical model single-particle states 
and/or occupation numbers of quasi-particle states as shown in ([97| and ((98|) . For 
simpler application and to increase computational speeds, these can be parameter- 
ized according to (jlOip . Spectroscopic factors for transitions from excited target 
states or the doubly averaged factors required in (|96p remain an open problem. 

The code SMARAGD (see Sec. I5.4.3P will also include a global DC treatment 
using an ADC model and energy-dependent spectroscopic factors, making use of a 
comb ination of t he above approaches intended to also yield consistency with the 
]j[p]y[[28 3 | 234 | 94 | g^f.]^ ADC approach aims at providing robust predictions de- 
spite of considerable differences between microscopic predictions.!^ Preliminary re- 
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Fig. 33. Comparison between averaged direct capture with energy-dependent spectroscopic factor 
(Sec. 15.53)1 . modified Hauser-Feshbach capture (Sec.[5..3.3ll, and standard Hauser-Feshbach capture 
(Sec. [53^ on even Sn isotopes (preliminary results) ."224, 

suits for this ADC are shown in Fig.|33l along with results from the HFM (Sec. l5.3."2l) 
and a modified HFM (Sec. I5.3.3|) . The final rate (or cross section) is the sum of 
the modified HFM value and the ADC one. Interestingly, for the isotopes shown 
here (except for = 92) this sum is approximated by the unmodified HFM result 
within a factor of 10. This is in accordance with Ref. 12841 (see figure 3 therein), 
where the DC contribution also almost replaces the standard HFM values. This 
shows that it seems justified to use unmodified HFM rates as crude estimate of the 
total rates for exotic nuclei. 

6. Conclusion 

There is no fast highway to improvements in predictions of reaction rates, not even 
a wide road. Historically, the fields of standard nuclear physics and nuclear data for 
applications have taken another direction. Therefore - instead of following a beaten 
track - rather a new, narrow path has to be driven step by step through a jungle of 
complications and details. This is only possible in a concerted effort of theory and 
experiment, the latter involving both large-scale rare-isotope production sites and 
smaller facilities. There is not one "most important" nucleus or "most important" 
reaction in Nuclear Astrophysics. Which nuclei, nuclear properties, and reactions 
are at the center of attention depends on the astrophysical process studied. There- 
fore, systematical studies are needed as well as information on specific nuclei and 
reactions. On the theory side large-scale studies of general trends and dependences 
are required as well as detailed predictions of individual nuclear properties. 

As shown above, the calculation of astrophysical rates, even when experimental 
information is present, involves a number of specialities not encountered in usual 
nuclear physics investigations. Thus, Nuclear Astrophysics is heavily relying on 
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advances in nuclear theory and experiment but also requires its own special devel- 
opments in theory and experiment which justify its existence as a separate field. 
Niiclcar physicists working at the boundary to astrophysics have to be aware of 
these special requirements and it is one of the aims of this work to having outlined 
a number of them. A further noteable fact is that it is necessary to not only point 
out special effects or important possible improvements but to actually apply them 
across the nuclear chart and produce large-scale sets of reaction rates which can 
be readily used by astrophysicists. This implies that they are made accessible in a 
form suitable to be implemented in astrophysical models. 

It should not be forgotten that the other essential aspect of Nuclear Astro- 
physics is the astrophysical modeling using reaction networks. Only in conjunction 
with this part of the field progress can be made. The models set the stage and define 
the ranges of conditions within which nuclear processes occur. Nevertheless, it can 
be treacherous to rely too strongly on a certain model. Reliable nuclear models and 
astrophysical reaction rates should cover a large range of possibilities and provide 
a sound base for pinpointing the sites of nucleosynthesis processes or even for dis- 
covering new types of nucleosynthesis in different astrophysical models. Moreover, 
the modifications of the cross sections and rates in a stellar plasma are an interest- 
ing topic in itself and warrant an independent study even without connection to a 
specific astrophysical site. 

The prediction of astrophysical reaction rates takes nuclear physics in a new 
direction and tests nuclear theory at the limits. Due to the finite number of nuclei, 
however, this is a finite task. This is also why parameterizations or phenomenolog- 
ical models may still have their justification, if designed in an appropriate manner. 
We do not have to extrapolate to infinity but within rather limited ranges of nu- 
cleon nimibers. It is reassuring that the overall abundance distributions obtained 
by combining several postulated nucleosynthesis processes are already closely re- 
sembling what we find in nature, even when explosive events and highly unstable 
nuclei are involved. This tells us that the most important properties seem to be 
described acceptably well. There are remarkable exceptions, however, both in ex- 
plaining abundance distributions (e.g., of the p- nuclei, the light s-process elements, 
isotopic anomalies in meteorites, the heaviest nuclei at the endpoint of the r-process, 
and many more) and in assigning astrophysical sites (e.g., to the r-process, proba- 
bly also to parts of the p-process). A detailed understanding and reliable reaction 
rates are also essential for using nuclear cosmochronometry to determine astronom- 
ical timescales. To go beyond previous estimates and reach a new level of detail, 
however, requires a large, dedicated effort of both experiment and theory. 
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